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ABSTRACT 


A  previous  report,  designated  Part  1,  with  the  same 
main  title,  presented  information  of  a  practical  nature,  with 
little  or  no  proofs,  derivations  of  equations,  or  explanations, 
except  that  sources  of  information  not  originated  by  the 
author  were  cited,  if  previously  published.  This  report. 
Part  2,  presents  the  proofs,  derivations,  and  explanations 
omitted  in  Part  1,  together  with  some  results  which  were 
either  previously  unpublished  or  are  now  out  of  print.  The 
basic  topics  covered  are  probability  of  detection  of  a  radar 
signal,  calculation  of  the  radar  pattern- propagation  factor, 
atmospheric  losses,  antenna  noise  temperature,  and  the 
relationship  of  range- prediction  errors  to  the  individual 
errors  in  the  range- equation  factors. 


PROBLEM  STATUS 

This  is  an  interim  report  on  a  continuing  NRL  Problem. 


AUTHORIZATION 

NRL  Problem  R02-55 
Project  RF05- 151-402-4011 


Manuscript  submitted  October  22,  1969. 


A  GUIDE  TO  BASIC  PULSE-RADAR  MAXIMUM-RANGE  CALCULATION 

PART  2  -  DERIVATIONS  OF  EQUATIONS,  BASES  OF  GRAPHS, 

AND  ADDITIONAL  EXPLANATIONS 


1.  INTRODUCTION 

Part  1  has  presented,  with  a  minimum  of  explanation,  the  equations,  definitions,  and 
data  required  for  calculating  the  maximum  range  of  a  pulse  radar.  In  Part  2,  the  deriva¬ 
tions  of  these  equations  will  be  given.  The  equations  and  methods  used  for  plotting  vari¬ 
ous  curves  presented  in  Part  1  will  also  be  given,  and  additional  explantions  of  some 
matters  will  be  furnished* 


2.  DERIVATION  OF  THE  RADAR  TRANSMISSION  EQUATION 

The  derivation  of  Eq.  (1),  Part  1,  the  radar  transmission  equation,  is  given  by 
Kerr  (1).  The  derivation  is  repeated  here  in  very  brief  fashion,  because  of  its  funda¬ 
mental  importance  to  the  subject  of  this  report. 

It  is  first  supposed  that  a  radar  transmitting  antenna  radiates  isotropically  —  that  is, 
with  an  equal  intensity  in  all  directions  of  space.  Then,  over  the  surface  of  an  imaginary 
transparent  sphere  surrounding  the  antenna,  which  is  located  at  the  center  of  the  sphere, 
assuming  a  lossless  propagation  medium,  the  power  density  (power  flowing  through  a  unit 
area  of  the  spherical  surface)  is,  by  invoking  the  principle  of  conservation  of  energy, 


where  Pt  is  the  total  power  radiated  and  R  is  the  radius  of  the  sphere.  (Its  total  surface 
area  is  4?rR2.) 

If  a  target  at  range  R  from  the  transmitting  antenna  intercepts  an  amount  of  power 
contained  in  an  area  a-  square  meters,  and  reradiates  it  isotropically,  the  reradiated 
power  density  returned  to  the  radar  antenna  will  be 


or,  substituting  for  St  the  right-hand  side  of  Eq.  (1), 


(4w)2  R4 


*Part  1  is  NRL  Report  6930,  dated  Dec.  23,  1969.  Its  subtitle  is  "Equations,  Definitions,  and  Aids 
to  Calculation," 
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An  actual  target  does  not,  of  course,  reradiate  isotropically;  but  the  standard  defini¬ 
tion*  of  radar  cross  section  is  based  on  Eq.  (2);  that  is, 

4ttR2  S_ 


In  other  words,  if  a  target  at  range  R  is  illuminated  with  a  power  density  St  watts  per 
square  meter  and  it  returns  to  the  radar  antenna  a  power  density  S„  its  radar  cross 
section  is  given  by  Eq.  (3).  (Free- space  propagation  conditions  are  here  assumed.)  Thus 
the  "radar  cross  section"  is  an  " equivalent  isotropically -reradiating  cross  section." 

If  the  radar  transmitting  antenna  has  a  power  gain  Gt,  in  the  direction  of  the  target, 
St  is  increased  by  the  factor  Gt;  that  is,  Eq.  (2a)  becomes 


S.  = 


ptGta 


r  (477)2  R4 

The  receiving  "capture  area"  of  an  antenna  of  gain  Gr  is 

G_  \2 


Ac  "  477 


(5) 


(Thus,  the  capture  area  of  an  isotropic  antenna  (G=  l)  is  A2/4tt.)  The  power  delivered 
to  a  receiver  by  a  lossless  receiving  antenna  of  capture  area  Ac  is,  by  definition  of  Ac, 


Pr  =  Sr  Ac  . 


(6) 


Making  the  substitutions  of  Eqs.  (4)  and  (5)  in  Eq.  (6)  gives 


Pt  Gt  Gr  a  X.2 
(4tt)3  R4 


(7) 


This  is  the  radar  transmission  equation  for  a  radar  and  target  located  in  free  space.  In 
the  actual  earth  environment,  St  and  Sr  are  further  modified  by  various  propagation  ef¬ 
fects,  such  as  refraction,  diffraction  (shadowing),  multipath  interference,  and  absorption. 
If  Sj.  and  S'  are  the  values  observed  in  the  actual  environment,  propagation  factors 
can  be  defined  such  that 


Ft2  =  s;/st 

(8a) 

Fr2  =  s;/sr  . 

(8b) 

(Actually,  a  propagation  factor  F  is  basically  defined  in  terms  of  the  ratio  of  the  electric 
field  intensity  observed  in  the  actual  environment  to  that  in  empty  space;  therefore,  F2 
applies  to  the  ratios  of  the  power  densities.)  The  resulting  radar  transmission  equation 
(expressed  in  terms  of  the  ratio  Pr/Pt)  is 


*The  word  "definition"  Is  in  italics  to  emphasize  the  fact  that  the  physical  unreality  of  the  isotropic- 
reradiation  assumption  is  not  objectionable,  since  definitions  are  to  a  large  extent  arbitrary.  The 
only  requirement  on  them  is  that  all  subsequent  uses  of  defined  terms  must  be  consistent  with  the 
definitions.  Additional  italics  in  the  paragraph  have  similar  motivation. 
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Pr  Gt  Gr  <r  \2  Ft2  Fr2 
(4t7)3  R4 


(9) 


This  completes  the  derivation.  (The  equation  is  given  by  Kerr  for  a  monostatic  radar 
using  the  same  antenna  for  transmitting  and  receiving.)  However,  thus  far  Gt  and  Gr  have 
been  assumed  to  be  antenna  gains  in  the  direction  of  the  target.  If,  following  Kerr,  they 
are  defined  as  gains  in  the  beam-maximum  directions,  then  the  reduction  of  gain,  if  any, 
in  the  target  direction  is  accounted  for  by  additional  factors  in  Ft  and  Fr.  These  factors 
are  then  called  pattern- propagation  factors.  It  has  been  shown  in  Part  1  how  this  basic 
equation  can  be  elaborated  to  obtain  radar  maximum  range  equations  in  terms  of  more 
readily  known  quantities,  and  for  more  complicated  situations. 


3.  SIGNAL- TO- NOISE  RATIO,  PROBABILITY  OF  DETECTION, 
AND  RELATED  SUBJECTS 


Single- Pulse  Detection 

The  curves  of  Figs.  4  through  7  of  Part  1  were  machine  computed  by  the  author  (and 
machine  plotted),  using  mathematical  results  derived  by  others  (with  minor  exceptions).* 


In  Part  1  it  was  explained  that  the  probability  of  detection  is  formally  given  by  the 
expression 

PsnO)  dv  -  (10) 

vt 

where  v  is  the  receiver  output  voltage  (after  all  predecision  processing,  including  inte¬ 
gration  or  averaging  of  successive  pulses  in  the  case  of  pulse  radar,  doppler  filtering  in 
the  case  of  cw  or  pulse-doppler  radar,  etc.)  The  quantity  vt  is  a  fixed  threshold  level  oi 
voltage  that  depends  (as  will  be  described)  on  the  allowable  rate  of  "false  alarms"  (deci¬ 
sion  that  a  signal  is  present  when  in  fact  there  is  no  signal).  The  function  psn  is  the 
probability  density  of  v,  which  is  a  random  variable  because  it  is  the  resultant  receiver 
output  when  both  signal  (s)  and  noise  (n)  are  present  in  the  input. 


The  probability  of  a  false  alarm  is 


(ID 


where  pn(v)  is  the  probability  density  of  the  detector  output  voltage  when  the  input  is 
noise,  with  no  signal  present.  The  function  psn  must  of  course  reduce  to  pn  when  the 
signal  power  has  zero  value. 

The  value  of  Pd,  Eq.  (10),  obviously  depends  on  the  value  selected  for  vt;  this  value 
is  chosen,  as  is  evident  from  Eq.  (11),  by  first  deciding  what  value  of  Pfa  is  desired,  or 
required,  and  then  solving  Eq.  (11)  for  the  value  of  vt  that  results  in  this  value  of  Pfa. 

Of  course,  this  solution  can  be  obtained  analytically  only  if  pn(v)  is  (in  some  sense)  an 
integrable  function. 


Figure  4  gives  values  of  Pd  as  a  function  of  the  signal- to- noise  power  ratio  S/N,  for 
detection  of  a  single  pulse  (no  integration),  and  a  linear-rectifier  detector  (sometimes 


*In  this  section,  reference  to  figures  is  to  those  of  Part  1  (NRL  Report  6930). 


4 


L.  V,  BLAKE 


referred  to  as  an  envelope  detector).  The  equations  for  pn(v)  and  psn(v)  for  this  case, 
as  has  been  shown  by  North  (2)  and  others,  are 


Pn<v> 


V  2 

p-vV 

n-  2 


_L  e-(v2/2cr2  +  S)  j 


,(£/=) 


where  a2  is  the  noise  power  at  the  detector  input,  and  S  is  the  signal-to-noise  power 
ratio  at  the  same  point.  These  equations  assume  a  Gaussian  distribution  for  the  ampli¬ 
tude  of  the  RF  noise  voltage.  This  assumption  is  valid  for  the  noise  sources  ordinarily 
present.  (l0  denotes  the  modified  Bessel  function  of  the  order  zero,  which  is  related  to 
the  ordinary  zero-order  Bessel  function  J0  by:  l0(x)  =  J0(ix);  i  =  ^FT.  Equation  (12) 
is  the  probability  density  function  of  a  Rayleigh  distribution,  and  that  of  Eq.  (13)  is  some 
times  referred  to  as  the  Rician  distribution  (after  S.O.  Rice  (3),  who  described  its  prop¬ 
erties  in  detail,  although  it  had  also  been  described  previously  by  North  and  others.) 

The  function  of  Eq.  (12)  is  integrable.  Therefore,  Eq.  (11)  can  be  solved  for  vt . 

The  result  is 

p  -  -t2/2-2  (141 


Therefore, 


-=  V- 2  l°ge  Pf£ 


The  ratio  vt/o-  is  the  ratio  of  the  threshold  voltage  to  the  predetection  rms  noise 
voltage.  It  is  more  practical  to  express  vt  in  relation  to  the  average  noise  voltage  at 
the  detector  output  (since  this  quantity  is  easily  measured;  it  is  the  dc  output  voltage  of 
the  detector).  Since 


Jro  _ _ 

vpn(v)  dv  =  <Jyj-n/2 


it  follows  that 


ut  =  —  =  v  4  l°gP  P f,A  ■ 


Robertson  (4)  has  normalized  the  threshold  voltage  differently;  the  relationship  of 
his  threshold  voltage,  ur,  to  the  one  used  here,  ut,  is 


V4  /n  -  1 


where  M  is  the  number  of  pulses  integrated  (m  =  1  for  Eq.  (17)). 
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Results  for  ut  and  also  ur,  obtained  by  digital  computer  solution  of  Eqs.  (17)  and 
(18),  are  listed  in  Table  1.  The  values  of  ur  in  Table  1  agree  with  those  reported  by 
Robertson  (4),  as  do  also  the  values  of  Pd(S)  plotted  in  Fig.  4.  This  agreement  is  note¬ 
worthy  because  the  methods  of  computation  were  entirely  different;  Robertson  employed 
a  series  representation  of  the  solution  for  Pd  ,  rather  than  direct  numerical  integration. 
His  results  are  presented  for  a  continuous  range  of  values  of  Pfa,  but  the  range  is 
smaller  than  the  one  of  Fig.  4.  Robertson's  curves  are  for  discrete  values  of  s  in  deci¬ 
bels  ranging  from  +6  to  +14. 


Table  1 

Values  of  Single-Pulse  Normalized 
Threshold  Voltage  for  a  Set  of  False- 
Alar  m  Probabilities 


Pfa 

ut 

Ur 

10"1 

1.7122 

1.3625 

lO-2 

2.4215 

2.7193 

10“3 

2.9657 

3.7605 

10“4 

3.4245 

4.6381 

10“5 

3.8287 

5.4114 

10“  6 

4.1941 

6.1105 

10“  7 

4.5301 

6.7534 

10“8 

4.8429 

7.3517 

10“9 

5.1367 

7.9137 

10-10 

5.4146 

8.4453 

10“ 11 

5.6788 

8.9509 

10“12 

5.9313 

9.4340 

10“13 

6.1735 

9.8973 

10-14 

6.4066 

10.3431 

10“15 

6.6315 

10.7733 

10-16 

6.8489 

11.1893 

These  values  of  ut  were  used  in  computing  the  curves  of  Fig.  4,  Part  1,  by  numeri¬ 
cal  integration  of  Eq.  (13),  with  v  replaced  by  u  =  v/v;  that  is, 


Pd(S) 


e-[(77U2/4)  +  S] 


IQ  (u  y/rr S)  du  . 


(19) 


The  numerical  integration  was  carried  out  to  seven- significant- figure  accuracy,  and  the 
results  were  machine- plotted.  The  original  plot,  of  which  Fig.  4  is  a  photographic  re¬ 
duction,  was  15  by  15  inches,  and  the  plotter  accuracy  is  specified  by  the  manufacturer 
to  be  ±0.005  inch  or  better. 

The  probability  scale  of  Fig.  4  resembles  a  Gaussian  probability  scale,  but  it  is  not; 
it  is  a  1/4-power-law  scale  which  has  similar  properties  but  is  more  convenient  for  ma¬ 
chine  computation.  This  scale  is  defined  by  the  relationships 


X 


>1/4 


-  P 


1/4 


,0.5 


1/4 


-  P 


1/4 


Pd  <  0.5  , 


(20a) 


and 
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2  - 


(1-Pd) 


1/4 


-  P. 


1  /  4 


0.5J 


-  P, 


1/4 


P,.  >0.5 


where  xro  is  the  midpoint  of  the  scale,  x  is  the  actual  coordinate  corresponding  to  a  prob¬ 
ability  Pd,  and  Pmin  is  the  minimum  value  of  Pd  (corresponding  to  x  =  0). 


It  is  interesting  to  note  that  the  Pd  curves  are  asymptotic  at  their  low  ends  to  the 
vertical  line  Pd  =  Pfa.  This  results  from  the  fact  that  by  definition  Pd  cannot  be  smaller 
than  Pfa;  in  other  words, 


lim  Pd(S)  =  Pfa  . 
s-*o 


(21) 


(This  result  is  indicated  by  the  fact  that  Eq.  (13)  becomes  identical  to  Eq.  (12)  when  S  =  0, 
because  io(0)  =  1.)  Of  course,  the  meaning  of  "probability  of  detection"  is  purely  formal 
in  this  region;  it  is  not  a  useful  concept  when  Pd  =  Pfa,  and  applications  of  Fig.  4  (and 
similar  curves  for  other  situations)  should  be  restricted  to  the  condition  Pd  »  Pfa,  How¬ 
ever,  the  fact  that  the  curves  of  Fig.  4  do  in  fact  become  asymptotic  to  the  Pd  =  Pfa  line 
is  a  check  on  the  correctness  of  the  computations. 


Detection  of  Two  Integrated  Pulses 

As  has  just  been  shown,  the  functions  pn(v)  and  psn(v)  are  exactly  known  for  detec¬ 
tion  of  a  single  pulse.  As  will  later  be  shown,  when  a  sufficient  number  of  pulses  are  in¬ 
tegrated,  approximate  formulas  which  are  very  accurate  can  be  given  for  pn  and  psn. 

For  two  or  just  a  few  pulses  integrated,  neither  the  exact  functions  nor  good  analytic  ap¬ 
proximations  can  be  found.  The  series  approximation  method  used  by  Robertson  is  prob¬ 
ably  the  best  method  of  solution.  However,  another  method  can  be  used  for  the  two-pulse 
case,  although  it  becomes  too  cumbersome  for  more  than  two  pulses.  This  method  was 
used  to  compute  the  two-pulse  value  used  in  Fig.  5,  Part  1.  It  was  used  partly  because 
it  is  a  method  that,  by  its  nature,  provides  an  interesting  check  on  more  analytical  meth¬ 
ods,  and  partly  because  the  author  wanted  to  see  if  it  could  be  done.  It  is  a  Monte  Carlo 
method,  utilizing  the  random-number  generator  which  is  an  on-line  function  of  the  NRL 
CDC-3800  computer. 

The  procedure  used  was  as  follows.  First,  the  function  psn,  expressed  in  the  form 
of  the  integrand  of  Eq.  (19),  was  "digitized"  by  dividing  the  values  of  v  into  increments 
of  0.01,  and  for  each  increment  the  value  of  psn  was  expressed  as  an  integer  equal  to  the 
rounded-off  value  of  100  psn  =  N.  Then  a  10 ,000- element  one-dimensional  array  was 
set  up  in  which  the  values  of  v  were  the  elements;  each  value  of  v  was  entered  into  the 
array  N  times,  where  N  is  related  to  psn(v)  as  stated  above.  Because  of  the  normaliza¬ 
tion  property  of  psn,  i.e., 


|  Psn(v)  dv  =  l  , 

the  division  of  the  range  of  v  into  increments  of  0.01,  and  the  multiplication  of  psn  by 
100,  the  resulting  10,000-element  array  contains  all  the  significant  part  of  the  psn 
function. 

An  array  of  this  type  was  set  up,  in  the  digital  computer,  for  each  value  of  S,  in 
decibels,  ranging  from  3.5  dB  to  15.5  dB,  in  0.5  dB  steps.  These  arrays  were  stored  on 
magnetic  tape. 
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To  compute  Pd(S),  using  these  arrays  and  the  Monte  Carlo  method,  the  array  for  a 
specified  value  of  s  was  buffered  into  the  computer  core  memory,  and  the  random  num¬ 
ber  function  was  used  to  select  an  integer  in  the  range  1  to  10,000;  each  such  integer  has 
equal  probability  of  being  selected.*  This  integer  is  then  used  as  the  array  index  to  se¬ 
lect  a  value  of  v,  designated  Vj.  This  represents  the  detected  voltage  for  the  first  echo 
pulse  of  a  two-pulse  pair.  The  procedure  is  then  repeated  to  obtain  a  second-pulse  volt¬ 
age,  v2.  The  sum  of  these  voltages,  vs  =  vt  +  v2,  then  represents  the  voltage  output  of 
a  two- pulse  postdetection  integrator. 

This  value  vs  was  then  compared  with  the  threshold  voltage,  vt,  corresponding  to  a 
desired  value  of  false-alarm  probability,  Pfa.  If  vs  >  vt,  a  detection  was  considered 
to  have  occurred;  if  vs  <  vt,  detection  did  not  occur.  This  "experiment"  was  repeated 
many  times,  and  a  running  account  was  kept  of  the  ratio  of  successes  to  the  total  number 
of  trials.  This  ratio,  when  taken  for  a  sufficient  number  of  trials,  is  a  good  experimental 
value  for  Pd(S).  The  trials  were  continued  until  this  experimental  value  stabilized,  as 
determined  by  the  following  criterion.  The  "running  average"  ratio  of  successful  detec¬ 
tions  to  total  trials  was  stored  at  the  end  of  each  sequence  of  1000  trials.  When  this 
average  for  ten  successive  times  (a  total  interval  of  10,000  "observations"  of  "pulse- 
pairs")  did  not  deviate  from  the  first  of  these  ten  values  by  more  than  0.00025,  the  run 
was  terminated,  and  the  value  of  Pd  at  the  termination  was  recorded.  Typical  numbers 
of  trials  required  for  satisfaction  of  this  criterion  ranged  from  about  30,000  to  150,000 
(a  "trial"  here  means  a  single  "pulse-pair"  addition).  The  0.00025  criterion  seems  to 
have  been  adequately  small  for  the  accuracy  needed.  It  was  noticed,  however,  that 
smaller  numbers  of  trials  were  needed  to  satisfy  the  criterion  at  very  high  values  of  Pd 
than  for  values  below  0.9.  After  the  fact,  it  seems  probable  that  the  accuracy  criterion 
should  be  made  more  stringent  at  very  high  (or  very  low)  values  of  Pd ,  but  the  value 
used  was  apparently  adequate  for  the  range  of  Pd  considered. 

About  20  minutes  of  computer  time  (on  the  CDC-3800)  were  used.  However,  it  was 
realized,  again  after  the  fact,  that  this  time  could  have  been  cut  considerably  if  a  test  for 
all  five  values  of  threshold  voltage  had  been  made  for  each  added  pair  of  randomly  se¬ 
lected  voltages  from  the  probability-density  arrays.  (Instead,  the  complete  experiment 
was  reperformed  for  each  value  of  threshold  voltage,  corresponding  to  a  specified  false- 
alarm  probability.) 

As  the  foregoing  discussion  indicates,  the  procedure  requires  determination  of  the 
threshold  voltage,  vt,  for  a  given  false-alarm  probability.  This  in  turn  requires,  in 
principle,  that  the  function  pn  be  given  analytically.  (It  is  not  required  that  the  function 
be  analytically  integrable,  as  it  happens  to  be  for  the  single-pulse  case.  If  pn  is  not 
analytically  integrable,  Pfa  can  be  found  by  numerical  integration,  for  a  specified  value 
of  vt.  However,  finding  vt  for  a  specified  Pfa  then  requires  either  interpolation  or 
iteration.) 

The  function  pn  for  the  two-pulse  case  cannot  be  obtained  as  an  analytic  function. 
Neither  is  it  feasible  to  find  Pfa  by  the  same  Monte  Carlo  method  that  has  just  been  de¬ 
scribed  for  finding  Pd .  This  method  is  feasible  only  for  probabilities  that  are  neither 
too  small  nor  too  large,  because  if  the  ratio  of  successes  or  failures  to  total  trials  in 
the  Monte  Carlo  experiment  is  too  small,  an  excessive  number  of  trials  is  required  for 
the  value  to  stabilize  —  that  is,  to  reduce  the  error  to  an  acceptable  value.  It  was  suc¬ 
cessfully  used  for  Pd  =  0.95,  but  it  would  probably  not  be  feasible,  for  example  at 


*The  author  is  aware  of  limitations  on  the  "randomness"  of  most  so-called  random-number  genera¬ 
tors,  but  analysis  of  the  method  and  of  the  results  indicates  that  these  limitations  did  not  affect  the 
validity  of  the  procedure  described  here. 
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Pd  =  0.999,  or  0.001.  For  Pfa  values  such  as  10"6  or  smaller,  it  would  be  out  of  the 
question.* 

However,  the  following  method  proved  to  be  practical.  The  probability  density  pn(2)> 
for  two  pulses,  is  related  formally  to  pn(l>  by  the  following  expression: 


p 


n( 


Pn(l)(x)  Pn(  1 )(  z  ~  x)  dx  ' 


(22) 


The  proof  of  this  formula  follows  directly  from  the  assumption  that  the  voltages  of  two 
successive  samples  of  the  detector  output  are  statistically  independent,  and  the  two  fol¬ 
lowing  propositions  of  probability  theory:  (1)  the  probability  of  the  simultaneous  occur¬ 
rence  of  two  independent  events  is  the  product  of  their  separate  probabilities,  and  (2)  the 
probability  of  any  single  event  is  equal  to  the  sum  of  the  probabilities  of  the  exhaustive 
and  mutually  exclusive  possibilities  of  its  occurrence.  The  product  of  the  two  terms  of 
the  integral  represents  the  application  of  proposition  (1),  and  the  integral  itself  is  an  ex¬ 
pression  of  proposition  (2). 

Since  pn(1)  is  given  by  Eq.  (12),  the  form  of  Eq.  (22)  in  the  present  instance,  taking 
cr  =  1  for  convenience,  is 


pn( 2)(z) 


=  J*  x  e_x2//2  (z-x)  e_(z_x)2,/2dx  . 


(23) 


After  some  manipulation  this  can  be  expressed  in  the  formt 

„z/2 


Pn(2)<»  =  e'z2/4 


(zV2-l)J  e'u2  du 


+  £  e-*2/2 
2 


(24) 


The  false-alarm  probability  is  now  related  to  this  result  by  Eq.  (11).  After  further 
manipulation,  this  yields: 


■  ; 


Pn(2)(Z>  dz 


=  0.5  e  Z‘  /2  +  J  e'z2/4  (z2/2-  1)  J 


z  /  2 


e-u  du 


d  z 


*Dr.  G.  V.  Trunk  of  the  NRL  Radar  Division  Analysis  Staff,  who  reviewed  portions  of  this  report, 
has  advised  the  author  that  he  has  used  the  Monte  Carlo  method  to  calculate  probability  of  detection 
for  n  pulses,  with  N  >  2  .  He  employed  the  characteristic  function  of  the  single-sample  noise  prob¬ 
ability  density  to  find  the  n -integrated-sample  density  function,  and  from  this  he  determined  the 
threshold  voltage  for  specified  false-alarm  probability.  He  then  employed  the  Monte  Carlo  method, 
as  described  here  for  two  pulses,  to  compute  probability  of  detection.  (This  method  is  basically 
applicable  to  any  number  of  pulses,  once  the  required  threshold  voltage  has  been  determined.) 
Trunk  has  also  found  it  feasible  to  use  the  Monte  Carlo  method  for  detection  probabilities  as  low 
as  0.001. 

tl  am  indebted  to  Dr.  G.  V,  Trunk  of  the  Radar  Division  Radar  Analysis  Staff  for  discovering  an 
algebraic  error  in  my  initial  derivation  of  this  result,  and  for  helpful  discussions  of  other  matters 
relating  to  the  computation  of  detection  probabilities. 
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The  numerical  solution  of  this  formula  requires  a  double  numerical  integration,  but 
this  is  easily  accomplished  by  the  digital  computer.  The  infinite  upper  limit  of  the  inte¬ 
gral  is  sufficiently  well  approximated  by  integrating  to  the  numerical  upper  limit  25. 

(The  result  obtained  for  the  overall  integral  of  Eq.  (25)  by  numerically  integrating  from 
0  to  25  differs  insignificantly  from  0.5,  which  is  the  analytically  correct  result,  since 
for  zt  =  0  in  Eq.  (25),  the  result  must  be  Pfa  =  1.) 

Since  Eq.  (25)  was  obtained  by  assuming  a  -  1,  the  average  voltage  output  for  one 
noise  sample  is  V2tT,  in  accordance  with  the  well-known  theorem  of  probability  theory  for 
the  average  of  the  sum  of  two  independent  random  variables.  Therefore,  the  threshold 
voltage  of  Eq.  (25)  is  related  to  threshold  ut  (normalized  to  the  average  voltage)  by  the 
relation 


Ut  =  Zt/yf2rr  . 


(26) 


Since  the  expression  derived  finds  Pfa  when  zt  is  specified,  and  actually  it  is  Pfa 
that  is  specified  in  probability  of  detection  calculations,  zt  was  found  for  specified  values 
of  Pf  by  an  interpolation  procedure.  The  results  for  a  range  of  values  of  Pfa  are  given 
in  Table  2,  in  terms  of  the  normalized  quantity  ut.  Also  the  equivalent  threshold  values 
u r ,  normalized  according  to  Robertson's  method  (4)  (Eq.  (18)  with  M  =  2),  are  given. 

Table  2 

Values  of  Normalized  Threshold  Voltage 
for  Two- Pulse  Detection 


Pfa 

ut 

ur 

1(T4 

2.644 

4.448 

10'6 

3.168 

5.866 

IO'8 

3.613 

7.069 

10-i° 

4.006 

8.133 

io- 12 

4.363 

9.099 

The  resulting  values  of  signal- to- noise  ratio  for  specified  probabilities  of  detection 
are  given  in  Table  3.  These  results  agree  well  (as  closely  as  reading  curves  can  affirm) 
with  Robertson's  results  within  the  range  of  values  of  Pfa  for  which  they  overlap,  and 
the  values  of  ur  given  in  Table  2  also  agree  well  with  those  published  by  Robertson. 

This  is  an  excellent  confirmation  of  both  calculations,  since  the  methods  were  radically 
different. 


Table  3 

Signal- to- Noise  Ratios  for  Detection  With  Specified 
Probability,  Two  Pulses  Integrated 


False- Alarm 
Probability 

Signal-to-Noise  Ratio  (dB)  for  Various 

Probabilities  of  Detection 

0.1 

0.25 

0.5 

0.75 

0.90 

0.95 

10' 4 

3.85 

5.48 

6.95 

8.22 

9.20 

9.74 

10" 6 

6.28 

7.51 

8.69 

9.71 

10.56 

11.02 

IO"  8 

7.87 

8.87 

9.88 

10.78 

11.52 

11.94 

10-io 

9.05 

9.95 

10.84 

11.65 

12.31 

12.68 

io- 12 

9.97 

10.79 

11.59 

12.33 

12.96 

13.31 
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In  the  Monte  Carlo  procedure  for  determining  Pd,  these  values  of  ut  were  rounded 
off  to  two  decimal  places,  since  that  was  the  limit  of  accuracy  set  by  the  number  of  ele¬ 
ments  in  the  digitized  arrays  of  detector  output  voltages. 


Detection  with  Many  Pulses  Integrated 

Although  the  method  that  has  been  described  for  detection  of  two  pulses  cannot  prac¬ 
tically  be  applied  for  even  three  pulses,  and  certainly  not  for  many  pulses,  certain  ap¬ 
proximations  can  be  made  for  the  functions  pn  and  psn  of  Eqs.  (10)  and  (11)  in  the 
many-pulse  case  which  become  very  good  if  the  number  of  pulses  is  sufficiently  large. 
These  approximations  are  described  by  North  (2).  They  make  use  of  the  central  limit 
theorem  of  probability  theory,  which  states  that  the  distribution  of  the  sum  of  n  inde¬ 
pendent  random  variables  tends  to  a  Gaussian  (or  normal)  form  as  n  tends  to  infinity, 
no  matter  what  the  distributions  are  for  the  individual  variables,  subject  to  some  fairly 
mild  conditions  that  are  well  satisfied  in  the  present  situation.  The  variances  and  aver¬ 
age  values  of  these  sum  distributions  are  found  by  applying  the  following  two  well-known 
theorems:  (1)  the  standard  deviation  crM  for  the  sum  of  M  independent  random  variables 
is  given  by 


°m2  =  £  ax2  -  (27) 

i=  1 

where  cri  is  the  standard  deviation  of  the  ith  variable;  (2)  the  average  of  the  sum,  xM,  is 
given  by 


xm  =  E  xi  •  (28) 

i=  1 

where  xi  is  the  average  value  of  the  ith  variable. 

The  desired  Gaussian  density  functions  for  the  detector  output  voltage  v  are  there¬ 
fore 


P(v) 


1  V>2/2cro2 

y/ln  ao 


(29) 


where  aD  and  v  are  found  by  applying  Eqs.  (27)  and  (28).  Since  a  steady  signal  and  con¬ 
stant  noise  power  are  assumed,  all  the  Vi's  and  c^'s  in  these  equations  are  equal;  hence 
for  M  pulses  integrated  the  results  are 


vM  =  M  Vj  (30) 

and 

°M  =  VM  cr1  .  (31) 

The  quantity  a0  in  Eq.  (29)  is  the  standard  deviation  of  the  detector  output  voltage, 
which  is  not  the  same  as  that  of  the  input  voltage,  denoted  o-  in  Eqs.  (12)  and  (13).  The 
relationships  of  v2  and  <rl  to  the  detector-input  quantity  a  for  the  noise- only  case  are 
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In 


r\J  77-/2 


and 


o-ln  =  a  \J  2  -  17/2  =  vlnV(4/7T)  -  1 

For  the  signal- and- noise  case  the  relationships  are* 


Vlsn  =  °V”  •  b  =  bvm  • 


CTlsn  =  vln  V(4/W)(l  +  S)  -  b2  , 


and 


b  =  e-s/2  [(1+  S)  I0(S/2)  +  S  •  IjCS/2)] 

(For  S  =  0,  b  =  1,  since  l0(0)  =  1  and  Ix(0)  =  0.) 

The  false- alarm  probability  is  then  given  by 


f  a 


V  2rr  < 


j-r 

CTMn  ^ 


■<V-VMn)  /2CTMn 


dv  , 


and  the  probability  of  detection  is 


?d  = 


-<v'vMsn>2/2crMSn 


dv 


Msn  vt 


(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


where  vt  is  the  selected  threshold  voltage. 

Integrals  of  the  type  shown  can  be  written  in  terms  of  the  error  function, 

u 

erf(u)  =— — I  e't2  dt  ,  (39) 

which  is  a  tabulated  function.  For  computer  programs,  standard  subroutines  for  erf  (u) 
are  available;  alternatively,  direct  numerical  integration  can  be  employed. 

In  terms  of  this  function,  and  with  quantities  in  the  arguments  of  the  error  function 
normalized  to  the  average  noise  output  (for  M  added  samples),  and  with  the  normalized 
threshold  voltage  written  as  u(t),  the  equations  for  false-alarm  probability  and  prob¬ 
ability  of  detection  become 


pfa  =  --  -  erf 
8  2  2 


ut  -  1 


V2(4/77-  1)/M 


(40) 


*See  S.  O.  Rice,  Ref.  3,  Eqs.  4.2-3  and  4.2-6,  pages  119-120.  The  notation  in(x)  in  Eq.  (36)  stands 
for  the  modified  Bessel  function  of  the  first  kind,  of  order  n.  For  nonnegative  integral  values  of  n, 
it  is  related  to  the  ordinary  Bessel  function  Jn(x)  by  the  formula  in(x)  =  i'n  jn  (ix);  i  =  \  -  1  • 
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This  equation  can  be  solved  for  ut  by  iteration,  or  by  interpolation  using  a  table  of 
the  error  function.  For  some  standard  values  of  Pfa  the  iteration  can  be  done  "in  gen¬ 
eral”  to  find  a  constant  kpfa,  such  that 

Ut  =  1  +  kpfa  V2(4/rr-  1)/M  .  (41) 


The  formula  for  kpfa  is 


kpfa  =  erf-»(l—  2  Pfa)  ,  (42) 

where  the  notation  erf-1  denotes  the  "inverse  error  function."  Some  useful  values  of 
kpfa,  which  were  obtained  by  interpolation  using  a  15-place  table  of  the  error  function,* 
are  given  in  Table  4. 


Table  4 

Values  of  Inverse  Error  Function,  kpfa, 
for  Some  Standard  Value  of  False  Alarm 
Probability, 


Pfa 

kpf  a 

10-4 

2.630 

10-6 

3.361 

10' 8 

3.968 

10-io 

4.498 

10-12 

• 

4.974 

The  corresponding  equation  for  Pd  is 


V2  [4(1  +  S)/t r  -  b2]/M  J 
Again,  a  constant  can  be  found  for  a  specified  value  of  Pd: 

kpd  -  erf-1  (1*2  Pj)  .  '  (44) 

In  terms  of  this  constant  Eq.  (43)  becomes 


-  pa 

V  2  [4(1+  S)/tt  -  b2]/M 

However,  this  equation  cannot  be  solved  explicitly  for  S,  because  b  is  a  function  of 
s.  If  the  explicit  expression  for  b  is  used  (Eq.  (36)),  it  becomes  evident  that  an  algebraic 
solution  is  not  possible.  Consequently,  iteration  is  required  for  numerical  solution  of 
Eq.  (45)  for  s,  just  as  it  would  be  if  Eq.  (43)  were  used  directly.  Nevertheless,  the  pro¬ 
cedure  of  finding  kpd  and  writing  Eq.  (45)  is  advantageous  when  solutions  for  several  or 


*For  small  values  of  pfa,  a  many-plaee  table  is  required,  such  as  the  15-place  "Tables  of  Probability 
Functions,"  VoL  1,  WPA  Tables,  U.S.  Govt.  Printing  Office,  1941.  These  tables  are  adequate  for 

Pf,  =  10-12. 
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many  values  of  N  are  desired  for  the  same  value  of  P^,  because  the  iterative  solution  of 
Eq.  (45)  does  not  involve  numerical  integration,  as  does  that  of  Eq.  (43).  Values  of  kpd 
for  some  standard  values  of  Pd  are  given  in  Table  5.  These  were  obtained  by  iteration 
using  the  CDC-3800  computer  and  an  error-function  subroutine.  It  is  noteworthy  that 

kpd(pd>  =  -kpd(l'pd)  -  <46) 

so  that  one  iteration  procedure  gives  values  of  kpd  for  two  values  of  Pd.  Also,  for 
Pd  =  0.5,  kpd  =  0  exactly. 


Table  5 

Values  of  Inverse  Error  Function,  kpd, 
for  Some  Standard  Values  of  Probability 
of  Detection,  Pd 


pd 

kPd 

0.01 

1.644976 

0.05 

1.163087 

0.10 

0.9061940 

0.25 

0.4769363 

0.50 

0.0000000 

0.75 

-0.4769363 

0.90 

-0.9061940 

0.95 

-1.163087 

0.99 

-1.644976 

This  "many-pulse"  method  for  computing  Pd  is,  of  course,  useful  only  if  some  idea 
can  be  given  of  the  errors  incurred  in  using  it  for  various  numbers  of  pulses.  Table  6 
compares  the  results  obtained  using  this  method  and  the  results  published  by  Robertson 
(4),  for  some  representative  cases.  Robertson's  results  were  calculated  for  numbers  of 
pulses  that  are  integral  powers  of  two  —  i.e.,  for  M  =  1,  2,  4,  8,  . . . ,  to  a  maximum  of 
8192  (=  213).  Table  6  makes  the  comparison  for  several  of  these  values  of  M,  for  two 
values  of  probability  of  detection,  and  for  those  values  of  false-alarm  probability  for 
which  the  two  sets  of  results  overlap.  It  was  difficult  to  read  Robertson's  curves  to  bet¬ 
ter  than  about  0.1  dB;  the  figures  given,  therefore,  sometimes  read  to  the  nearest  0.1  dB, 
sometimes  to  one  more  decimal  place.  The  decibel  discrepancy  was,  however,  in  all 
cases  rounded  off  to  the  nearest  0.1  dB.  It  is  evident  that  the  agreement  improves,  as 
would  be  expected,  as  the  number  of  pulses  increases  —  because  the  assumption  of 
Gaussian  density  curves  for  pn  and  psn  (Eqs.  (10)  and  (11))  is  more  nearly  correct  as 
the  number  of  pulses  integrated  increases  (in  accordance  with  the  central  limit  theorem). 
It  is  probably  a  good  rule  of  thumb  that  the  method  is  virtually  exact  for  numbers  of 
pulses  in  the  thousands,  very  good  for  numbers  in  the  hundreds  (error  of  order  0.1  dB), 
and  in  error  by  roughly  0.5  dB  at  about  ten  pulses  integrated. 

It  will  be  evident  to  the  reader  that  Robertson's  results  have  been  used  as  a  check  on 
the  correctness  of  results  obtained  for  each  of  the  three  computational  procedures  em¬ 
ployed  (for  one,  two,  and  many  pulses  integrated).  It  might  thus  seem  that  it  would  have 
been  simpler  in  the  first  place  to  plot  the  curves  desired  using  Robertson's  results,  or 
to  use  his  method  for  computing  the  plotted  points.  As  a  matter  of  fact,  the  procedures 
described  here  were  developed  before  Robertson's  paper  was  published.  Moreover,  his 
paper  does  not  give  the  details  of  his  computational  procedure.  Quite  apart  from  these 
practical  factors,  however,  there  are  several  justifications  for  carrying  out  the  proce¬ 
dure  that  has  been  described.  First,  just  as  Robertson's  results  confirm  the  correctness 
of  the  results  obtained  here,  so  also  do  these  results  confirm  Robertson's.  Second,  the 
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Table  6 

Comparison  of  North  and  Robertson  Methods  for 
Several  Values  of  Parameters 


No.  of 
Pulses 

False- Alarm 
Probability 

Probability 
of  Detection 

Signal-to-Noise 

Ratio  (dB) 

Discrepancy 

(dB) 

North’s 

Method 

Robertson's 

Curves 

8 

10'4 

0.5 

2.14 

2.70 

0.9 

4.26 

4.60 

10‘6 

0.5 

3.42 

4.1 

0.9 

5.23 

5.8 

0.6 

32 

10’4 

0.5 

-1.25 

-0.95 

0.9 

0.54 

0.70 

10'6 

0.5 

-0.08 

0.28 

0.4 

0.9 

1.44 

1.72 

0.3 

128 

10'4 

0.5 

-4.45 

-4.3 

0.2 

0.9 

-2.88 

-2.8 

0.1 

10‘6 

0.5 

-3.33 

-3.19 

0.1 

0.9 

-2.02 

-1.90 

0.1 

4096 

10-4 

0.5 

-12.13 

-12.1 

0 

0.9 

-10.80 

-10.8 

0 

lO-6 

0.5 

-11.06 

-11.02 

0 

0.9 

-9.96 

-9.95 

0 

10'8 

0.5 

-10.33 

-10.1 

0.2 

0.9 

-9.38 

-9.45 

0.1 

curves  that  have.been  plotted  here  cover  a  range  of  values  of  the  false-alarm  probability 
not  available  in  Robertson's  published  paper.  Finally,  the  methods  employed  here  are 
inherently  interesting,  for  various  reasons;  they  illuminate  some  of  the  theory  of  detec¬ 
tion—  especially  the  Monte  Carlo  procedure  for  the  two-pulse  case,  since  it  is  based 
partly  on  a  mathematical  analysis  and  partly  on  a  statistical  experiment.  It  is  inter¬ 
mediate  between  a  purely  mathematical  solution  and  an  experimental  one  in  which  an 
actual  threshold  device  with  real  signals  and  noise  would  be  employed. 


Curve  Plotting  Procedure 

The  curves  of  Figs.  5  through  22,  Part  1,  were  plotted  by  machine.  The  plotter  used 
for  these  curves  was  the  Gerber  Model  875,  which  plots  to  an  accuracy  of  ±0.005  inch  or 
better,  with  overall  plot  dimensions  up  to  5  by  8  feet.  The  plots  were  made  to  a  size  of 
16  by  20  inches,  and  photographically  reduced  for  use  in  this  report.  The  coordinate 
grid  was  also  done  by  the  machine,  to  avoid  any  problems  of  registration  which  would 
arise  if  the  plotting  were  done  on  ready-made  coordinate  paper.  (The  coordinate-axis 
numbering  and  labeling  was  also  done  by  the  machine.) 

Since  points  of  the  curves  were  computed  for  a  limited  set  of  values  of  the  number  of 
pulses,  a  plotting  subroutine  was  used  which  computes  interpolated  points  to  produce  a 
smooth  curve,  in  much  the  same  manner  that  a  draftsman  would  use  a  French  curve. 

(The  subroutine,  written  by  the  author,  is  named  FCURVE.) 
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To  avoid  possible  inaccuracy,  the  "many-pulse”  method  was  used  for  Figs.  5  through 
10  only  for  very  large  numbers  of  pulses.  After  some  experimenting  it  was  found  that 
curves  of  best  accuracy  (as  judged  by  comparison  with  Robertson's  curves)  could  be 
produced  by  computing  five  points,  at  M  =  1,  2,  8000,  9000,  and  10,000  pulses.  The  latter 
three  points  establish  a  slope  at  the  many-pulse  end  of  the  curve,  and  the  one-  and  two- 
pulse  points  establish  a  slope  at  the  low  end.  The  FCURVE  subroutine  provides  the  re¬ 
quired  transition  between  these  slopes.  This  method  produced  a  better  agreement  with 
Robertson's  results,  in  the  region  of  overlap,  than  did  computation  of  additional  points 
with  the  many-pulse  method.  The  agreement  finally  produced  is  as  good,  for  all  numbers 
of  pulses,  as  the  readability  of  the  curves. 

In  accordance  with  theory,  the  slope  of  these  curves  is  asymptotic  to  -1.5  dB  per 
octave  in  the  small- signal  region  (S/N  «  l),  and  to  -3  dB  per  octave  in  the  large- signal 
region  (S/N  »  l).  That  is,  S/N  <=  M-1  for  S/N  »  l  and  S/N  M_1/2  for  S/N  «  1.  Rob¬ 
ertson  has  noted  that  his  results  are  similarly  in  agreement  with  these  asumptotic  values. 


Results  for  Fluctuating  Signals 

Figures  6  and  7  of  Part  1,  which  are  for  the  Swerling  Cases  1  and  3  types  of  fluc¬ 
tuating  signals,  were  plotted  in  the  manner  just  described,  with  points  computed  using  a 
computer  program  due  to  Fehlner  et  al.  of  the  Johns  Hopkins  University  Applied  Physics 
Laboratory,  and  described  in  Ref.  5.  It  is  essentially  a  confirmation,  with  some  correc¬ 
tions,  of  the  results  of  Marcum  (6)  and  Swerling  (7).  It  was  modified  by  the  author  and  by 
Stanley  Gontarek*  to  permit  it  to  run  on  the  CDC-3800  computer  instead  of  on  an  IBM 
machine  for  which  it  was  written.  It  was  further  modified  by  the  author  to  compute  re¬ 
sults  with  false-alarm  probability  as  an  input  parameter,  rather  than  Marcum's  false- 
alarm  number.  Figures  6  and  7  are  also  plotted  in  the  same  format  as  Fig.  5,  with 
number  of  pulses  integrated  and  signal- to- noise  ratio  as  coordinates  and  false-alarm 
probability  as  a  family  parameter.  (The  curves  published  by  Robertson  have,  as  coor¬ 
dinates,  probability  of  detection  and  false-alarm  probability;  the  family  parameter  is 
signal- to- noise  ratio.  Fehlner's  curves  have  signal- to- noise  ratio  and  probability  of 
detection  as  coordinate  axes,  with  Marcum's  false-alarm  number  as  the  family  param¬ 
eter.) 

In  the  literature  of  radar  range  calculation,  there  is  sometimes  defined  a  "fluctuating- 
signal  loss,"  which  is  applied  when  the  nonfluctuating- signal  visibility  factor  has  been 
used  in  a  range  calculation  and  it  is  desired  to  convert  the  result  to  apply  to  the  fluctuating- 
signal  case.  As  a  casual  inspection  of  Figs.  5  through  7  indicates,  this  concept  is  not 
tenable.  The  "loss"  is  sometimes  a  gain,  and  it  cannot  be  related  to  any  one  of  the  other 
parameters.  The  relationship  of  the  signal- to- noise  ratio  required  for  detection  in  the 
steady- signal  case  to  that  for  fluctuating- signals  does  not  permit  a  simple  correction  to 
be  made.  The  correct  procedure  is  to  calculate  the  range  using  a  visibility  factor  ob¬ 
tained  from  the  appropriate  curve  of  Figs.  5  through  7,  or  some  similar  curve  appro¬ 
priate  to  the  applicable  fluctuation  model. 


Calculation  of  Radar  Range  for  Specified  Target- Cross- 
Section  Measurement  Accuracy 

In  Part  1,  radar  "maximum  range"  equations  are  given  for  specified  values  of  prob¬ 
ability  of  detection  and  specified  false-alarm  probability.  The  foregoing  sections  of 
Part  2  have  detailed  the  theory  of  this  maximum- range  concept,  which  is  the  appropriate 


*Naval  Air  Systems  Command. 


16 


L.  V.  BLAKE 


one  for  most  radar  applications.  However,  in  certain  uses  of  radar,  it  is  desired  to  know 
the  range  at  which  the  signal-to-noise  ratio  will  reach  a  specified  value,  rather  than  the 
range  for  a  specified  probability  of  detection.  These  are  primarily  measurement  uses  — 
that  is,  application  in  which  the  radar  is  being  used  to  measure  some  characteristic  of 
the  radar  target,  such  as  its  radar  cross  section.  In  these  applications,  the  appropriate 
equation  is  Eq.  (6)  of  Part  1,  with  the  "max"  and  "min"  subscripts  removed. 

The  question  to  be  answered  in  these  applications  is:  how  does  the  accuracy  of  the 
measurement  depend  on  the  signal-to-noise  ratio? 

Two  cases  (at  least)  can  be  considered.  The  first  is  that  in  which  the  desired  meas¬ 
urement  is  to  be  made  on  the  basis  of  a  single  received  pulse  (or  a  single  sample  of  a 
CW  signal,  the  sample  time  being  equal  to  or  less  than  the  reciprocal  of  the  receiver 
noise  bandwidth).  In  the  second  case,  the  measurement  is  made  after  averaging  a  num¬ 
ber  of  pulses  or  samples  of  the  signal  and  noise.  For  the  first  case,  assuming  also  that 
only  reasonably  small  errors  are  to  be  permitted  (so  that  only  large  signal-to-noise 
ratios  are  considered),  the  following  analysis  applies.* 

If  the  receiving  system  were  completely  noise  free,  the  measured  receiver  output 
signal  would  be  an  accurate  measure  of  the  signal  power  at  the  receiver  input.  The 
presence  of  noise,  however,  causes  the  receiver  output  voltage  at  a  specific  instant  to 
be  randomly  somewhat  greater  or  less  than  the  value  which  would  be  observed  under 
noise-free  conditions.  Accordingly,  a  measurement  of  this  voltage  made  at  any  particular 
instant  (i.e.,  for  any  single  received  pulse)  will  be  in  error  by  some  amount  that  cannot 
be  predicted.  However,  the  statistics  of  this  error  can  be  calculated,  as  a  function  of  the 
input  signal-to-noise  power  ratio.  The  basic  theory  from  which  such  calculation  may  be 
made  has  been  developed  by  Rice  (3)  and  others.  A  linear-rectifier  detector  is  assumed. 
Also  assumed  is  a  linear  receiver-amplification  characteristic.  Although  in  practice 
there  may  be  some  nonlinearity,  correction  for  this  nonlinearity  can  be  made  in  process¬ 
ing  the  data. 

The  following  notation  will  be  used: 

P  =  amplitude  of  sine-wave  RF  input  signal  voltage, 

R  =  rms  value  of  input  RF  noise  voltage, 

S  =  P2/2R2,  input  signal-to-noise  power  ratio, 
v  =  receiver  output  voltage  (video  pulse  voltage). 

Rice  has  shown  (Ref.  3,  Eq.  (4.2-5),  p.  120)  that,  assuming  a  linear-rectifier  detec¬ 
tor,  for  large  signal-to-noise  ratios  the  average  output  voltage  is  asymptotically  given  by 

v  =  C  (P  +  R2/2P  +  R4/8P3+  ...)  (47) 

and  that  the  mean- square  value  (Rice,  Eq.  (4.2-6))  is 

v2  =  C2  (P2  +  2R2)  .  (48) 


j|5 

This  analysis  is  basically  taken  from  the  author’s  Appendix  A  of  Ref.  8,  but  has  been  extended  to 
apply  to  target-cross-section  measurement. 
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(In  these  expressions  c  is  a  constant  proportional  to  the  voltage  gain  of  the  receiver.) 
Therefore  the  variance  a2  of  v  is  (neglecting  terms  of  order  higher  than  R2/2P2) 

a2  =  ^2  -  v2  =  C2R2(1-  R2/2P2)  ,  (49) 

corresponding  to  Rice's  Eq.  (4.2-7).  Hence  the  standard  deviation  is 

o-  3  CR(1  -  R2/2P2)1/2  .  (50) 

Moreover,  Rice  has  shown  that  for  sufficiently  large  signal- to- noise  ratios,  the  dis¬ 
tribution  of  v  is  approximately  normal,  meaning  that  the  probability  density  is  given  very 
nearly  by 


p(v)=—^ —  exp  [-  (  v  -  v  )  2/  2cr2  ]  .  (51) 

V2  tt  a 

The  probability  that  the  observed  voltage  deviates  from  the  average  value,  v,  by 
more  than  a  specified  amount  Av  is  therefore  given  by 

?  =  Prob.  ( |  v  -  v  |  >  Av)  =2  j'  P(v)  dv  .  (52) 

v  +  Av 

For  a  given  deviation,  say  Av  =  ka,  the  fractional  error  of  the  voltage  measurement, 
Av/v,  from  Eqs.  (47)  and  (50),  is  (again  neglecting  higher-order  terms)  approximately 

Av  _  ko;  ^  kP  (  1  ~  R2/2P2)  1/2 
V  V  (P  +  R2/2P+ R4/8P3) 

-  kR  (  1  -  R2/2P2) 1/2 
P  (1  +  R2/2P2  +  R4/8P4) 


1/2 

k  (1-1/4S) 

V2S  (1  +  1/4  S+  1/32S2) 


(53) 


in  which  S,  as  previously  noted,  is  the  input  signal-to-noise  power  ratio.  Inasmuch  as 
it  has  already  been  assumed,  in  deriving  this  expression,  that  s  is  much  larger  than 
unity,  the  term  in  brackets  will  ordinarily  be  approximately  equal  to  unity,  so  that  for 
most  calculations  the  simple  formula 


Av  k 
v  vTs 


(54) 


may  be  used. 

This  analysis  relates  the  accuracy  of  the  signal  voltage  measurement  (after  detec¬ 
tion)  to  the  predetection  signal-to-noise  ratio.  In  many  applications,  although  this  volt¬ 
age  (v)  is  the  directly  measured  quantity,  the  ultimate  objective  is  measurement  of  the 
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radar  cross  section*  of  the  target,  crt.  This  is  accomplished  through  the  following  rela¬ 
tion  (which  applies,  as  may  be  deduced  from  Eq.  (47),  when  the  detector  is  linear  and  the 
signal-to-noise  ratio  is  large,  because  P  is  proportional  to  at); 

crt  =  kcv2R4/F4  ,  (55) 

where  kc  is  a  calibration  constant,  usually  evaluated  by  observing  a  "standard  target" 
whose  radar  cross  section  is  known— for  example,  a  conducting  sphere;  r  is  the  radar 
range;  and  F  is  the  pattern- propagation  factor. 

From  Eq.  (55)  a  relationship  between  the  fractional  error  of  measurement  of  v  and 
the  corresponding  error  of  measurement  of  o-t  can  be  written.  The  procedure  is  first  to 
write  the  difference  equation: 


crt  +  A at  =  kc(v+ Av)2R4/F4  .  (56) 

where  A<rt  is  an  increment  of  <rt  corresponding  to  the  voltage  increment  Av.  This  leads 
to  the  result 


A<xt 


(57) 


Equation  (52)  can  be  rewritten  as  follows  in  terms  of  a  new  variable  t  =  (v-  v)/&: 


The  quantity 


9  = 


\JTtt 


e-t2/2  dt  =  1 


Av/cr 


t/277 


Av/cr 


s_t  '/2  dt 


E(k)  - 


e-t2/2  dt  , 


(58) 


(59) 


sometimes  known  as  the  error  integral,  is  a  tabulated  function,  and  is  also  readily  eval¬ 
uated  by  numerical  integration  using  a  computer.  The  latter  method  was  used  here.  For 
various  values  of  k  and  S,  9  was  found  from  Eq.  (58),  Av/vwas  found  from  Eq.  (54),  and 
A crt/at  was  found  from  Eq.  (57).  Figure  1  is  a  plot  of  9  as  a  function  of  Acrt/o-t  with  s 
(in  decibels)  as  a  parameter. 

For  Case  2,  in  which  many  pulses  are  averaged,  the  signal-to-noise  ratio  is  generally 
small,  typically  much  less  than  1.  Again  assuming  a  linear  detector,  the  average  detec¬ 
tor  output  when  there  is  no  signal  is  given  by  Eq.  (32),  and  with  a  signal  present  it  is 
given  by  Eq.  (34).  These  two  averages  are  related  by  the  factor  b,  Eq.  (36).  The  ques¬ 
tion  is  then,  what  is  the  limitation  on  the  accuracy  imposed  by  the  fluctuations  in  the 
measured  values  of  the  quantities  vsn  and  vn? 


*It  is  unfortunate  in  the  following  discussion  that  the  symbol  a  is  used  for  both  target  cross  section 
and  for  the  standard  deviation  of  signal  and  noise  voltages.  This  is  done  because  a  is  the  standard 
symbol  for  both  of  these  quantities.  The  subscript  t  will  be  used  whenever  a  target  cross  section 
is  meant. 
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Fig.  1  -  Probability  of  specified 
errors  in  single-pulse  measure¬ 
ment  of  radar  cross  section 


The  analysis  that  follows  is  based  on  the  assumption  that  the  signal-to-noise  power 
ratio,  S,  is  very  small  (much  less  than  1),  the  opposite  of  the  assumption  made  for  the 
one-pulse  measurement  case.  (The  intermediate  situation  is  much  more  difficult  to 
analyze,  but  the  results  of  an  analysis  would  necessarily  be  intermediate  between  those 
of  the  two  limiting  cases.)  Postdetection  integration  of  M  pulses  is  assumed. 

The  small- argument  approximation  for  the  exponential  and  modified  Bessel  functions, 
appropriate  when  S  «  i,  result  in  the  following  first-order  approximation  of  Eq.  (36): 

b  =  1  +  S/2  .  (60) 

Applying  this  result  to  Eq.  (34)  and  solving  for  s  gives 


(61) 


For  a  single  pulse,  Eqs.  (33)  and  (35)  relate  o-ln  and  crlsn  to  vln.  Applying  Eq  (60) 
to  Eq.  (35)  shows  that  for  S  «  1,  aln  **  crlsn,  so  that  only  Eq.  (33)  is  needed: 


1  sn 


rln  =  v In  V4/7T  -  1 


(62) 


If  M  pulses  are  integrated, 


and 


(See  Eqs.  (30)  and  (31).)  Therefore 


aMn  =  °-ln 


Mvj  . 


°Mn 


vMn 

Vm 


V4/77  -  1 


(63) 

(64) 


(65) 


and 
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2Dm 

5Mn 

Again  a  difference  equation  can  be  written: 


(66) 


S  +  AS  = 


2  (Dm  +  adm) 

^Mn 


(67) 


Here  vMn  is  treated  as  a  constant,  since  it  can  be.  measured  with  much  greater  accuracy, 
ordinarily,  than  can  vM,  or  %,  because  it  can  be  averaged  over  a  much  longer  time 
period. 

Dividing  by  s  on  the  left  and  by  2DM/vMn  on  the  right  gives 


AS  _  ADm 
S  '  Dm 


(68) 


Since  S  is  directly  proportional  to  the  target  cross  section  at ,  this  means  that  the  per¬ 
centage  error  in  measurement  of  aq  is  equal  to  the  corresponding  error  in  measurement 
of  When  the  number  of  pulses  integrated  is  large,  the  probability  density  function  of 
Dm  is  normal  (gaussian).  (This  is  approximately  true  for  as  few  as  ten  pulses  integrated.) 
Therefore  the  procedure  followed  in  the  one-pulse  error  analysis  (Eq.  (52))  can  be  ap¬ 
plied.  The  equation  corresponding  to  Eq.  (53)  (making  use  of  Eqs.  (65)  and  (66),  and  re¬ 
calling  that  the  standard  deviation  of  DM  is  o-Mn)  is 


A(Tt  ADm  kaMn  2k  - 

— —  =  -  =:  -  =  —  V4/7T  -  1 

D,,  l\,  S^M 


(69) 


Combining  all  of  the  numerical  factors  and  omitting  the  intermediate  expressions  gives 


Aoq 


cr 


t 


1.0454 


k 

S  VNT 


(70) 


Following  the  same  procedure  described  for  Fig.  1  results  in  a  similar  plot,  Fig.  2,  with 
Syff  as  the  parameter  instead  of  S. 


The  foregoing  analysis  applies  to  the  pulse-radar  measurement  of  a  point-target 
radar  cross  section.  Another  situation  will  now  be  treated  briefly  because  it  applies  to 
a  practical  problem  of  some  importance  —  the  measurement  of  the  signal  power  returned 
from  a  volume  distribution  of  independently  random  scattering  elements.  This  situation 
arises  for  example  in  measurement  of  the  Thomson  scatter  echo  from  the  ionosphere, 
which  is  used  to  determine  the  electron-density  profile.  Other  examples  are  echoes 
from  rain  and  from  "chaff"  (a  cloud  of  metallic  reflectors  used  as  a  radar  countermeas¬ 
ure),  The  problem  will  be  discussed  in  terms  of  the  Thomson  scatter  echo,  but  the  re¬ 
sults  are  applicable  to  the  other  examples  also.  This  echo  signal  is  well  known  to  be  a 
Gaussian-noise-like  signal,  so  that  its  effect  is  to  increase  the  level  of  the  receiver  noise 
already  present.  The  statistics  of  the  combined  signal  and  noise  are  therefore  identical 
to  those  of  the  receiver  noise,  and  Eqs.  (29),  (32),  and  (33)  apply.  Moreover,  if  the  Thom¬ 
son  scatter  signal  power  at  the  input  of  a  linear  detector  is  wTS,  the  standard  deviation 
is  ctts  =  \/wxs,  and  the  standard  deviation  of  the  total  power  at  the  same  point  is 
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Fig.  2  -  Probability  of  specified 
errors  in  many-pulse  measure¬ 
ment  of  radar  cross  section 
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cr 


sn 


+  CT 


2 

n 


(71) 


where  an  is  the  standard  deviation  of  the  receiver  noise.  Therefore, 


TS 


CT 


2 

n 


(72) 


From  Eq.  (32),  the  average  detector  output  voltages  vsn  and  vn  can  be  substituted  on 
the  right-hand  side,  giving 


"TS  =  T  KV  Vn2) 


(73) 


Thus  the  desired  measurement  of  wTS  can  be  made  by  measuring  the  observable  quanti¬ 
ties  Vsn  and  Vn,  the  average  detector  output  voltages  when  the  Thomson  scatter  signal  is 
present  and  absent. 

Again  it  will  be  assumed  that  vn  can  be  measured  much  more  accurately  than  vsn, 
so  that  the  error  in  determining  wTS  is  due  primarily  to  the  error  of  measurement  of 
Vsn.  The  relationship  of  the  fractional  error  in  measurement  of  wTS  to  that  of  vsn  can 
as  before  be  obtained  by  writing  a  difference  equation,  which  after  some  manipulation 
yields 


Awts 

WTS 


2vc„  Av  „  +  Av^ 

sn  sn  sn 


(74) 


Making  use  of  Eqs.  (32)  and  (71),  the  denominator  on  the  right  can  be  rewritten  in 
the  form 


o-jS/o-n2  is  the  input  signal- to- noise  power  ratio.  Eq.  (74)  thus  becomes 


in  which  S  = 
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AwTS  _  / S  +  l\ 
WTS  "  \  S  ) 


(76) 


This  equation  applies  for  "single  sample"  measurement  of  vsn»  In  this  case  the 
probability  density  function  is  not  gaussian  (Eq.  (29));  it  is  the  Rayleigh  function  (Eq.  (12)). 
The  Rayleigh  probability  density  function  for  a  single  sample  of  the  signal-plus- noise 
output  is 


P(v) 


JL 1  p-7rv2/4v2 

2v2 


(77) 


(This  equation  is  equivalent  to  Eq.  (12)  except  that  it  is  in  terms  of  v  rather  than  a, 
since  vsn  is  the  quantity  to  be  measured;  see  also  Eq.  (32).)  The  probability  of  making 
a  measurement  error  exceeding  Av  is  therefore 


9=  l- 


v- Av 


77V 

2v2 


e-7TV2/4v2 


dv  . 


(78) 


This  integration  is  performed  by  making  the  change  of  variable  u  =  7rv2/4v2,  where¬ 
upon 


where 


and 


9  =  1  - 


1  -  e 


-L„ 


-L. 


+  e 


(v-Av)2 


t  77  —  As  77 

!  =  W* (  1,>  =  7 


1  +  2 


Av 


(79) 


(80) 


(81) 


subject  to  the  constraint  Av  =  v  (since  the  lower  limit  in  Eq.  (78)  cannot  be  less  than 
zero).  When  numerical  results  using  this  formula  are  examined,  it  is  found  that  values 
of  Av/v  corresponding  to  Awts/wxs  of  the  order  of  unity  give  values  of  9  greater  than  0.5, 
regardless  of  how  large  s  becomes  (since  the  factor  (S  +  l  )/s  approaches  unity  asymp¬ 
totically  as  S  tends  to  infinity).  This  means  that  a  one-pulse  Thomson  scatter  measure¬ 
ment  is  subject  to  an  unacceptably  large  probable  error,  as  is  well  known  by  practitioners 
of  this  method  of  ionospheric  measurement. 


The  measurement  is  customarily  made,  in  Thomson  scatter  work,  by  averaging  many 
samples  which  are  statistically  independent  (corresponding  to  the  signal  returns  from 
many  successive  pulses  transmitted  by  the  radar).  The  distribution  of  this  average  is 
asymptotically  Gaussian,  by  the  central  limit  theorem  (discussed  preceding  Eq.  (27)). 
Therefore  f  is  calculated  from  the  formula  (similar  to  Eq,  (58)) 


\f2rt  J0 


9=1- 


e-t2/2  dt  , 


(82) 
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where  k  is  Av/o-.  The  relationship  between  k  and  Av/vsn  is  (applying  Eq.  (32) 

Av  kcr  _  k 

^sn  ^sn  V^72 


(83) 


If  M  statistically  independent  samples  are  averaged,  the  standard  deviation  of  the 
distribution  of  vsn  is  reduced  by  the  factor  compared  to  that  of  a  single  sample,  as 
indicated  by  Eq.  (31).  Consequently,  the  probability  decreases  for  exceeding  a  given 
value  of  Av/v  as  M  increases.  That  is,  since  the  standard  deviation  of  vsn  is  reduced  by 
the  factor  VM,  so  also  is  the  value  of  Avsn  corresponding  to  a  probability  9  reduced  by 
the  same  factor. 

Consequently  Eq.  (76)  can  be  rewritten  in  the  form 


(84) 


where  (Awts/wts)  '  denotes  the  fractional  error  in  measuring  wTS  when  M  samples  are 
averaged  and  Avsn/vsn  is  the  one-sample  value.  If 


1 

2  VM 


«  1  , 


(85) 


the  second  term  in  the  brackets  of  Eq.  (84)  can  be  neglected,  and  the  equation  can  then  be 
rewritten  in  the  form 


Vm 


(86) 


Therefore  if  9  is  calculated  for  the  one-sample  value  of  Avsn/vsn,  from  Eq.  (82), 
and  the  results  plotted  as  a  function  of  2  (Avsn/vsn),  the  curve  can  be  interpreted  as 
applicable  to  the  parameter  vOvT  [S/(S+  i>]  (AwTS/wTS) ,  as  shown  in  Fig.  3,  provided  that 
M  is  actually  appreciably  greater  than  1  (because  Eq.  (82)  does  not  really  apply  for 
M  =  1).  (The  prime  symbol  on  (Awts/wxs)  has  been  dropped,  but  is  implied,  in  the  label¬ 
ing  of  the  abscissa  scale.) 


After  (Awts/wts)  has  been  found  using  this  curve,  for  specified  values  of  S,  M,  and 
9,  a  test  should  be  made  to  insure  that  Eq.  (85)  is  satisfied.  The  equivalent  relation  in 
terms  of  (AwTS/wTS)  is 


(This  relation  can  always  be  satisfied  by  taking  m  sufficiently  large.) 

For  some  specified  value  of  9,  a  simple  equation  can  be  written  for  (Awts/wts).  For 
example,  if  9  =  0.001,  the  abscissa  of  the  curve  is  5.25.  Therefore,  for  9  =  0.001, 

Awts  _  5J_25  /S  +  l\ 
wts  VM"  \  S  / 


(88) 
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Fig,  3  -  Probability  of  specified 
errors  in  Thomson-seatter  signal- 
power  measurement 


Equation  (89)  is  satisfied  in  this  case  if  VM  »  1.3.* 

As  this  discussion  shows,  there  may  well  be  some  practical  applications  for  which 
values  of  M  are  used  that  do  not  satisfy  this  criterion.  In  such  cases  the  more  accurate 
expression,  Eq.  (84),  must  be  used,  but  then  no  simple  curve  relating  9  and  Awxs/wts  can 
be  drawn.  It  becomes  necessary  to  draw  a  family  of  curves  with  M  as  a  parameter,  or 
to  select  a  specific  value  of  9  and  for  this  value  draw  a  curve  for  [S/(S+  l)]  (AwTS/wTS) 
as  a  function  of  m  .  (Such  curves  would  be  analogous  to  those  of  Figs.  5  through  7,  Part  1.) 
Table  7  shows  how  the  values  of  [S/(S+  1 )]  (AwTS/wTS)  vary  with  M  for  the  particular  case 
9  =  0.001.  The  last  two  columns  show  how  the  result  of  this  exact  computation  compare 
with  the  approximate  result  that  would  be  obtained  from  Fig.  3.  That  is,  if  the  number  in 
column  3  is  called  x,  and  the  number  in  column  4  is  y, 

y  =  100  (1  -  5.  2514/x)  .  (89) 

As  indicated,  the  approximate  result  is  good  enough  for  most  purposes  for  M  >  1000.  The 
numbers  in  the  third  column  are  to  be  compared  with  the  value  5,2514  obtained  by  the 
approximate  method  for  this  value  of  9. 

Table  8  relating  9  and  k  is  useful  for  making  calculations  of  this  type.  The  values 
were  machine  computed  from  Eq.  (82)  for  the  specified  values  of  9  by  an  iterative 
method. 


Measurement  of  Minimum  Detectable  Signal 

Radar  receivers  are  often  tested  by  measuring  the  "minimum  detectable  signal,"  as 
a  check  on  their  level  of  performance.  Actually  a  noise  factor  or  noise  temperature 
measurement  is  the  best  measure  of  receiver  performance,  but  the  minimum -signal 
measurement  can  be  made  somewhat  more  easily,  and  is  thus  a  more  practical  method 


*It  is  to  be  noted  here  that  the  test  should  be  made  in  terms  of  y'M ,  ratherthan  in  terms  of  m,  to  avoid 
confusion  as  to  the  meaning  of  the  symbol  »,  For  example,  if  %/m  >>1.3  is  assumed  to  mean  the 
same  as  M  »  1 . 3  2 ,  then  if  » is  taken  to  mean  "two  orders  of  magnitude"  for  the  square-root  form 
of  the  relationship.  It  should  be  interpreted  as  "four  orders  of  magnitude"  when  applied  to  the  sec¬ 
ond  form. 
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Table  7 

Variation  of  [S/(S  +  l )]  (AwTS/wTS)  With  M  for  9  =  0.001 


M 

/  S  \  /Awts\ 

\S  +  lj  \  WTS  j 

Col.  2  x  VM 

Error  Using 
Fig.  3 

d) 

10 

2.3500 

7.431 

29 

100 

0.5941 

5.941 

12 

1,000 

0.1730 

5.469 

4.0 

10,000 

0.05320 

5.320 

1.3 

100,000 

0.01668 

5.273 

0.41 

1,000,000 

0.005258 

5.258 

0;13 

Table  8 

Relation  of  9  and  k 


9 

k 

0.001 

3.2908 

0.01 

2.5758 

0.1 

1.6449 

0.5 

0.67449 

of  routine  performance  testing.  In  fact,  if  the  measurement  is  made  properly,  it  is  in 
effect  an  indirect  measurement  of  the  noise  factor.  As  it  is  usually  carried  out,  however,- 
it  lacks  the  precision  usually  attained  in  a  properly  made  noise- factor  measurement. 

The  lack  of  precision  is  associated  with  the  subjective  nature  of  the  measurement  —  its 
dependence  on  the  judgment  and  experience  of  the  person  who  performs  it  and  makes  the 
decision  as  to  what  signal  level  is  "just  detectable.” 

Moreover,  considerable  confusion  sometimes  exists  as  to  the  value  of  minimum  de¬ 
tectable  signal  that  should  be  observed  if  the  receiver  is  operating  properly  (i.e.,  if  its 
noise  factor  is  as  good  as  it  should  be) .  The  purpose  of  the  following  discussion  is  to 
derive  the  correct  formula  for  computing  this  value.* 

As  can  be  deduced  from  Eqs.  (8)  and  (9)  of  Part  1,  the  following  equation  relates  the 
minimum  detectable  signal  power  (watts)  and  other  parameters  of  the  radar  system: 

Pr(min)  =  k  Ts  VQ  CB/T  .  (90) 

where  k  =  1.38  *  10"23  watt-sec/deg  is  Boltzmann’s  constant,  v0  is  the  visibility  factor, 
CB  is  the  bandwidth  correction  factor  (Eq.  (17)  and  Fig.  1,  Part  1),  and  r  is  the  radar 
pulse  length  in  seconds. 

The  minimum- signal  measurement  is  usually  made  with  a  signal  generator  connected 
to  the  receiver  input.  Therefore  the  system  noise  temperature  is 


*This  section  of  the  report  is  based  on  an  informal  paper  written  by  the  author  in  1962  in  the  Search 
Radar  Branch,  Radar  Division,  NRL. 
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Ts  =  Tg  +  Te  ,  (91) 

where  Tg  is  the  Kelvin  temperature  of  the  signal- generator  output  impedance,  and  Te  is 
the  receiver  noise  temperature.  Te  is  related  to  the  receiver  noise  factor  Fn  by  Eq. 
(41),  Part  1: 


Te=To(Fn-l)-  (92) 

where  T0  is  the  IEEE  reference  temperature  for  noise-factor  measurement,  290  °K.  If 
Tg  =  T0,  as  is  usually  true,  at  least  within  the  precision  of  other  factors  in  the  measure¬ 
ment,  Eq.  (91)  becomes 


T 


S 


F„  T„ 


(93) 


and  Eq.  (90)  becomes 


Pronin)  =  kFnT0V0CB/r  .  (94) 

It  is  customary  to  express  Pr  (min)  in  decibels  relative  to  1  milliwatt.  The  equation 
in  these  terms,  with  r  expressed  in  microseconds  and  with  the  numerical  values  of  k 
and  T0  inserted,  is 

Pr(min)(dBm)  =  +  V0(dB)  +  CB(dB)  +  Fn(dB)  "  10  %sec  ' 

For  most  well-designed  radar  systems,  CB  =  l,  or  CB(dB)  =  o.  The  usual  method 
of  making  the  measurement  is  to  connect  the  receiver  output  to  an  oscilloscope,  corre¬ 
sponding  to  a  radar  A- scope  display.  A  calibrated  pulse  signal  generator  is  connected 
to  the  receiver  input,  with  the  pulse  length  and  repetition  frequency  set  to  the  appropriate 
values.  The  visibility  factor  to  be  used  in  Eq.  (95)  is  then  given  by  a  curve  of  the  type  of 
Fig.  8,  Part  1.  As  mentioned  in  the  figure  caption,  this  curve  is  for  0.5  probability  of  de¬ 
tection  and  for  an  unspecified  false-alarm  probability.  A  human  observer  performing  a 
minimum-detectable- signal  measurement  seldom  if  ever  attempts  to  evaluate  these  pa¬ 
rameters.  This  is  another  reason  for  the  relative  imprecision  of  such  measurements, 
compared  to  noise-factor  measurement.  Experience  has  shown  that  the  typical  observer 
probably  evaluates  Pr  (min)  more  nearly  at  the  0.9  probability  level  than  at  the  0.5  level. 
The  0.9  probability  values  are  approximately  2  dB  higher  than  the  0.5  values,  for  pulse 
repetition  frequencies  below  about  1000.  (At  higher  repetition  frequencies  the  0.9  and  0.5 
probability  curves  gradually  come  together  and  tend  asymptotically  to  approximately  the 
constant  value  v0(dB)  = -8.) 

A  possible  source  of  error  in  the  measurement  is  the  use  of  too  low  a  sweep  speed 
on  the  oscilloscope.  As  discussed  in  Part  1,  a  loss  of  signal  detectability  occurs  if 
sr  <  d,  where  s  is  the  sweep  speed  in  millimeters  per  second,  r  is  the  signal  pulse 
length  in  seconds,  and  d  is  the  cathode- ray- tube  spot  diameter  in  millimeters.  There¬ 
fore,  care  should  be  taken  to  use  an  adequate  sweep  speed.  Usually  a  value  of  st  eaual 
to  or  greater  than  1  mm  is  sufficient. 

Sometimes  in  measurement  of  Pr(min)  the  receiver  is  left  connected  to  the  antenna, 
and  the  signal  from  the  signal  generator  is  fed  in  through  a  directional  coupler.  With 
this  arrangement,  the  system  noise  temperature  in  Eq.  (90)  is  given  by  the  following 
equation,  rather  than  by  Eq.  (9 1) : 


Ts  =  (Ta+Tr)/Lr  +  Te  . 


(96) 
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In  this  equation  Ta  is  given  by  Eq.  (37)  and  Tr  by  Eq.  (40)  of  Part  i.  Lr  is  the  receiving- 
transmission- line  loss  factor.  Equation  (96)  is  based  on  the  assumption  that  the  direc¬ 
tional  coupler  is  located  at  the  receiver  input.  The  attenuation  of  the  coupler  must  of 
course  be  considered  in  evaluating  the  signal  power  at  the  receiver,  but  it  does  not  affect 
the  noise  calculation  if  the  coupling  factor  (signal  attenuation)  is  large  —  e.g.,  20  dB  or 
greater.  Evaluation  of  Ta  and  Tr  is  discussed  in  Part  1  and  also  later  in  this  report. 
(Equation  (36)  of  Part  1  would  apply  instead  of  Eq.  (96)  if  the  directional  coupler  were 
located  at  the  antenna  terminals  instead  of  at  the  receiver  terminals.)  The  equation 
corresponding  to  Eq.  (95)  for  this  situation  is 

Pr(min)(dBm)  =  “138.6  +  VQ ( dB ^  +  CB(rfB^  +  10  log  Ts  -  10  log  T^sec  .  (97) 

4.  CALCULATION  OF  PROPAGATION  FACTOR 

In  Part  1,  a  method  of  calculating  the  pattern- propagation  factor  F  that  results  from 
sea- reflection  interference  when  earth  curvature  is  significant  is  described.  This  method 
is  taken  from  Ref.  1.  It  is  based  on  the  use  of  parameters  S  and  T,  from  which  the  fac¬ 
tors  J,  K,  and  D  are  calculated.  These  factors  are  graphed  in  Fig.  17  of  Part  1.  The 
graphs  were  machine  plotted  from  digital- computer  solution  of  the  appropriate  equations. 

In  this  section  the  method  of  calculating  these  graphs  will  be  described,  and  some  addi¬ 
tional  topics  relating  to  the  pattern-propagation  factor  will  be  discussed. 


Computer  Solutions  for  J,  K,  and  D 

The  basic  equations  for  computing  J,  K,  and  D  as  functions  of  s  and  T  are  given  in 
Ref.  1,  Sec.  2.13,  p.  113  ff.  However,  the  solutions  given  there  are  not  explicit.  There¬ 
fore  the  equations  on  which  the  computer  solutions  were  based,  resulting  in  the  graphs 
presented  in  Part  1,  will  be  presented  here  briefly.  The  starting  point  is  the  assumption 
that  the  parameters  s  and  T  are  given;  they  are  calculated  from  Eqs.  (57)  and  (58)  of 
Part  1. 


In  each  case,  the  calculation  is  done  through  intermediate  variables  Sj 
are  related  to  S  and  T  by  the  equations 

and  S2 ,  which 

S  =  (S1  +  S2T)/(1  +  T) 

(98) 

and 

S2(l  “  Si2)  -  SjT (  1  -  S22)  . 

(99) 

Eq.  (99)  can  be  solved  for  S2, 

giving 

S2  = 

\/(l  -  Sj2)2  +  4Sj2  T2  +  Sx2  -  1 

/2SJ  . 

(100) 

The  equations  for  J,  K,  and  D  in  terms  of  these  quantities  are 

J  =  (1- Sj2  )  (1- S22)  ,  (101) 

( 1  -  s  2 )  +  T2(l-  S  2) 

1/  _  _  _ * _ _  _ ± _ 


1  +  T2 


(102) 
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and 


D  = 


l  + 


4Sj2S2T 


-1-1/2 


S(l-S2)(l  +  T) 


(103) 


A  suitable  plotting  procedure  would  be  to  start  with  arbitrarily  chosen  values  of  St 
and  T,  calculate  S2  and  S  from  Eq.  (100)  and  (98),  then  find  J,  K,  and  D  from  Eqs.  (101), 
(102),  and  (103).  It  was  decided  however,  to  plot  contours  of  constant  J  and  K.  There¬ 
fore  the  procedure  actually  followed  for  J  was  to  select  arbitrary  values  of  J  and  T,  then 
to  iterate  Eq.  (101)  with  respect  to  s1}  with  the  appropriate  substitution  for  S2  given  by 
Eq.  (100).  In  this  way  the  required  values  of  S2  and  S2  (through  Eq.  (100)  were  found. 

Then  S  was  found  from  Eq.  (98).  A  similar  procedure  was  followed  for  K, 

A  plot  of  D,  however,  was  done  differently,  because  it  was  decided  that  a  more  read¬ 
able  graph  would  result  if  D  were  plotted  (on  a  logarithmic  scale)  as  a  function  of  s  with 
T  as  a  parameter.  Therefore  the  procedure  for  each  selected  value  of  T  was  to  go  through 
a  suitable  range  of  values  for  sx  to  allow  plotting  the  curve  of  D,  finding  the  values  of  S 
from  (in  order  of  their  use)  Eqs.  (100)  and  (98). 

The  plot  of  Fig.  18  of  Part  1  was  made  by  computing  J  and  K  as  just  described;  D 
was  also  computed  by  the  method  described  for  j  and  K  —  that  is,  by  the  iterative  method 
starting  with  the  appropriate  values  of  S  and  T.  Then  F  was  computed  from  (in  order  of 
their  use)  Eqs.  (59),  (60),  (56),  (62),  and  (61)  of  Part  1. 


Plots  of  v  and  u  for  the  Diffraction  Region 


The  quantities  V  and  U  plotted  in  Fig.  19  of  Part  1  are  based  on  Figs.  2.19  and  2,20 
of  Ref.  1.  The  curves  for  V  as  a  function  of  range  X  in  natural  units  were  computed  from 
the  basic  equation  (see  Ref.  1,  p.  122) 


VdB  =  20 


log 


10 


2  \/rrX  e"2  02  X 


(104) 


The  curves  for  u  were  simply  replotted  from  Fig.  2.20  of  Ref.  1,  since  the  computation 
is  quite  complicated;  it  is  given  by  Eq.  (428),  p.  109,  of  Ref.  1. 


The  "Fine  Lobe"  Case 

In  Part  1,  formulas  were  developed  for  calculating  the  pattern  propagation  factor  F. 
An  implicit  assumption  was  made  that  the  target  moves  sufficiently  slowly  in  range  to 
make  this  calculation  meaningful.  That  is,  it  was  assumed  that  if  the  radar  scans  in 
azimuth,  the  distance  that  the  target  moves  in  the  time  of  one  scan  is  such  that  F  does 
not  change  greatly.  In  terms  of  lobe  diagrams  such  as  Figs.  14  and  18  of  Part  1,  it  was 
assumed  that  the  range  change  ar  in  the  time  of  one  scan  is  a  small  fraction  of  the  range 
extent  of  a  lobe. 

If  this  assumption  is  not  correct,  then  the  calculation  of  F  as  described  is  not  mean¬ 
ingful.  If  a  plot  of  F  as  a  function  of  range  were  made  from  experimental  data  under  such 
circumstances,  the  result  would  not  be  a  plot  having  the  appearance  of  Fig.  18,  Part  1.  In 
effect,  F  behaves  in  this  case  like  a  randomly  fluctuating  quantity  (stochastic  variable). 
The  blip-scan  ratio  discussed  in  Part  1,  and  defined  there  by  Eq.  (86),  is  likewise  a 
stochastic  quantity  rather  than  an  analytic  variable.  The  meaningful  quantity  is  its 
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average  value  </<,  which  can  be  calculated,  as  shown  by  Alderson  (9).  The  following 
derivation  is  excerpted  from  Alderson's  report. 

The  general  equation  for  F,  for  the  sea- reflection-interference  case,  as  given  in 
Part  1,  Eq.  (61),  can  be  written  in  the  form 

F(<f)  =  |f  (0)Vl  +  X2  +  2x  cos  g\  ,  (105) 

where  f  (9)  is  the  antenna  pattern  factor,  x  is  the  magnitude  of  the  "effective  reflection 
coefficient,"  and  g  is  the  phase  difference  between  the  direct  and  reflected  rays.  As  a 
target  approaches  in  range  R,  g  changes  rapidly  while  9  and  x  change  relatively  slowly. 
Therefore,  4>  can  be  found  by  averaging  over  g.  The  angle  g  is  given  by 


f  = 


4vh  sin  9 
X 


J(S,T)  +  4>  , 


(106) 


in  which  J  is  the  correction  factor  discussed  earlier  (Eq.  101),  h  is  the  radar  antenna 
height,  X  is  the  wavelength,  0  is  the  target  elevation  angle,  and  is  the  phase  angle  of 
the  sea  reflection  coefficient.  (Since  many  of  the  equations  of  Part  1  are  in  terms  of 
"mixed"  systems  of  units,  it  should  be  noted  that  in  this  equation  h,  X,  and  R  must  all  be 
in  the  same  units  of  length.)  The  correction  factor  J  is  equal  to  unity  when  earth  curva¬ 
ture  can  be  neglected. 

The  nature  of  the  variation  of  < p  with  range  when  F  =  1  (free  space)  is  due  primarily 
to  the  target- cross- section  fluctuation,  which  has  been  discussed  in  Part  1.  The  statistics 
of  this  fluctuation  depend  on  the  nature  of  the  target.  In  the  absence  of  specific  target- 
fluctuation  information,  the  Rayleigh  probability- density  function  is  the  one  that  is  usually 
assumed  to  apply  to  the  signal  voltage  at  the  receiver.*  The  corresponding  density  func¬ 
tion  for  the  target  cross  section  is  given  by  Eq.  (88),  Part  1,  and  the  blip- scan  ratio  is 
then  given  by  Eq.  (89),  Part  1.  Therefore,  the  average  value  of  i/i,  averaged  over  the  full 
range  of  g  ( 2rr  radians)  is 


J27T 

(107) 

where  r50  is  the  range  at  which  \p  =  0.5  with  F  =  1  (free  space).  (This  average  is  equiv¬ 
alent  to  averaging  \p  over  one  full  lobe  of  the  interference  pattern.)  Alderson  numerically 


*The  Rayleigh  probability  density  function  for  a  voltage  v  is 

o  “  (  v/v  so  ) 
p(v)  =  2  In (  2  )  ( v/Vg 0 )  2 

where  v50  is  the  median  value  of  the  voltage.  The  more  usual  representation,  which  is  entirely 
equivalent,  is  given  by  Eq.  (77).  This  voltage  density  function  corresponds  to  a  target-cross¬ 
section  fluctuation  whose  density  function  is 


p(<0 


ln  2  -o-ArS0 

^  5  0 


This  is  a  negative-exponential  density  function  (not  the  Rayleigh  function);  but  a  target  whose 
cross-section  fluctuation  is  thus  characterized  is  commonly  called  a  "Rayleigh  target,"  because  of 
the  resulting  signal-voltage  fluctuation. 
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integrated  this  equation  to  plot  </»  as  a  function  of  r/r50,  for  various  values  of  the  param¬ 
eter  x,  Eq.  (105),  (Alderson  equated  this  parameter  to  D,  the  divergence  factor,  as¬ 
suming  that  p0,  r,  and  f  (52)/f  (01)~  see  Eq.  (62),  Part  1  —  were  all  equal  to  unity  as  is 
often  approximately  true.  Alderson's  curves  may  be  interpreted  for  the  more  general 
situation  by  simply  reading  x  in  place  of  his  D.)  The  resulting  curves  are  shown  as 
Fig.  1  of  Ref.  9.  They  show  that  for  D  =  0,  corresponding  to  free  space  (no  reflected 
wave),  the  variation  of  >A  with  range  R  is  steeper  than  for  nonzero  values  of  D.  In  other 
words,  the_ additional  signal  fluctuation  caused  by  a  "fine  lobe"  structure  results  in  lower 
values  of  at  ranges  appreciably  less  than  Rso,  and  larger  values  of  y5  at  ranges 
greater  than  Rs0,  compared  to  the  free-space  values.  See  Appendix  (p.  47). 

As  noted  in  Part  1,  these  results  are  calculated  without  taking  into  account  the  addi¬ 
tional  fluctuation  of  the  echo  signals  caused  by  the  presence  of  noise  in  the  receiving 
system;  the  more  exact  analysis,  including  the  effect  of  noise,  has  been  made  by  Mallett 
and  Brennan  (10).  However,  Alderson's  results  are  valid  when  the  minimum  detectable 
signal-to-noise  ratio  is  relatively  large  (true  when  detection  is  based  on  observation  of 
a  single  pulse  or  integration  of  only  a  few  pulses);  moreover,  the  general  conclusions 
drawn  from  his  curves  are  also  valid  in  the  more  general  case.* 

Alderson  (9)  has  applied  this  type  of  analysis  to  the  problem  of  analyzing  the  effects 
of  ship  pitch  and  roll  on  radar  performance ,  as  a  function  of  the  antenna  vertical  beam- 
with  in  relation  to  the  pitch  and  roll  angles. 


The  Effect  of  Diffuse  Forward  Scattering 

In  Part  1,  the  reduction  of  the  effective  reflection  coefficient  caused  by  roughness  of 
the  reflecting  surface  was  discussed,  and  a  formula  (Eq.  (56),  Part  1)  was  given  for  com¬ 
puting  the  ratio  p/p0,  where  p  is  the  magnitude  of  the  rough- surface  reflection  coeffi¬ 
cient,  and  p0  is  the  smooth- sea  coefficient. 

In  this  context,  p  may  be  called  the  specular  reflection  coefficient.  Of  the  total 
power  incident  upon  a  unit  area  of  the  sea  surface,  a  fraction  p2  is  specularly  reflected 
and  another  fraction  ( l  -  p2)  is  transmitted  through  the  surface  and  eventually  absorbed. 
The  remaining  fraction,  equal  to  p2  -  p2 ,  is  reflected  diffusely .  That  is,  there  is  no 
coherent  or  analytical  relationship  of  the  phases  of  this  portion  of  the  reflected  field  at 
points  in  space  separated  by  appreciable  distances.  Moreover,  in  the  case  of  the  sea 
surface,  this  phase  relationship  changes  randomly  with  time. 

The  diffusely  reflected  field  can  thus  be  regarded  as  a  randomly  varying  quantity  in 
both  space  and  time.  Its  existence,  however,  means  that  the  total  field  at  any  point  in 
space  is  really  the  vector  sum  of  three  components  —  the  direct- path  field,  the  specularly 
reflected  field,  and  the  diffuse  (randomly  varying)  field. 

Because  the  diffusely  reflected  field  can  be  regarded  as  the  vector  sum  of  reflections 
from  numerous  independently  related  reflecting  elements  of  the  sea  surface  (11),  the  am¬ 
plitude  distribution  of  this  field  can  be  assumed  to  be  Gaussian.  It  thus  has  a  noiselike 
character,  and  its  combination  at  the  receiver  input  with  the  direct-path  and  specularly- 
reflected  components  will  result,  at  the  output  of  a  linear  detector,  in  a  Rician  distribu¬ 
tion  of  signal  amplitude  (see  Eq.  (13)).  The  diffuse  forward- scattered  signal  therefore 


*  Alderson  notes  that  a  procedure  similar  to  his  was  described  by  T.  E.  Phipps  and  S.  L,  Waleszczak 
in  OEG  (Chief  of  Naval  Operations)  Study  473,  29  May  1952  (Confidential).  He  also  credits  A.  G. 
Ferris,  then  of  NRL,  with  the  conception  of  Fig.  1,  Ref.  9. 
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produces  signal  fluctuation,  in  addition  to  any  fluctuation  that  may  already  be  present 
from  causes  previously  discussed. 

The  random  field  due  to  diffuse  reflection  will  at  some  times  and  at  some  points  in 
space  add  to  the  field  resulting  from  interference  of  the  direct  and  specularly  reflected 
waves,  and  at  other  times  and  points  it  will  subtract.  On  the  average,  however,  since 
this  field  represents  additional  reflected  energy,  it  will  increase  the  field  strength,  and 
increase  the  probability  of  detection  of  a  radar  target.  In  the  lobe  maxima,  where  the 
specularly  reflected  signal  is  strongest,  the  relative  effect  of  the  diffuse  reflection  is 
least.  Its  effect  increases  as  the  ratio  p/p0  decreases  —  i.e.,  as  the  sea-wave  height  gets 
larger  in  relation  to  the  product  of  the  radar  wavelength  and  the  cosecant  of  the  grazing 
angle. 

However,  it  must  not  be  supposed  that  the  power  ratio  of  the  diffuse  to  the  specular 
field  at  a  radar  target  can  be  deduced  from  the  ratio  p/p0  —  that  is,  the  power  ratio  at 
the  target  is  not  p02/p2  -  1,  because,  although  this  does  represent  the  ratio  of  the  dif¬ 
fusely  to  specularly  reflected  power  per  unit  area  of  the  sea,  the  diffusely  reflected 
power  is  scattered  in  all  directions  while  the  specularly  reflected  power  is  concentrated 
in  the  forward  direction.  On  the  other  hand,  diffusely  reflected  signal  power  can  reach 
the  target  from  a  wide  area  of  the  sea,  if  the  radar  beam  is  wide,  while  the  specularly 
reflected  signal  is  simply  proportional  to  p2  regardless  of  the  radar  beamwidth  as  long 
as  it  is  not  narrower  (in  the  reflected- ray  direction)  than  the  angular  width  of  the  first 
Fresnel  zone.  Therefore  the  quantitative  relationship  of  the  diffuse  and  specular  signal 
components  is  quite  complicated. 

The  effect  of  the  diffuse  component  will  be  most  significant,  of  course,  in  the  "nulls" 
of  a  smooth- sea  interference  pattern.  Some  further  details  of  the  effect  are  given  by 
Blake  (11)  and  subsequently  by  Beckmann  and  Spizzichino  (12)  and  others.  The  effect  is 
mentioned  by  Kerr  (1),  page  535. 


5.  PROPAGATION- MEDIUM  LOSSES 
Tropospheric  Absorption 

In  Part  1,  curves  are  given  for  the  loss  due  to  absorption  of  radar  waves  in  the  nor¬ 
mal  troposphere.  Figure  21  gives  the  two-way  (radar)  absorption  for  total  traversal  of 
the  troposphere;  these  results  would  apply,  for  example,  for  radar  targets  in  outer  space, 
such  as  satellites  and  other  space  objects.  Since  ionospheric  absorption  is  not  included, 
these  curves  apply  fully  only  at  frequencies  at  which  ionospheric  absorption  is  negligible — 
generally,  above  about  1000  MHz  in  the  daytime,  and  above  perhaps  500  MHz  at  night.  At 
lower  frequencies  they  apply  for  targets  that  are  above  the  troposphere  but  below  the 
ionosphere. 

For  radars  used  to  detect  aircraft  and  other  targets  in  the  troposphere,  curves  are 
needed  which  give  the  absorption  as  a  function  of  range  and  angle.  Such  curves  are  pre¬ 
sented  in  Part  1  as  Fig.  22,  with  the  radar  frequency  as  a  parameter.  The  method  by 
which  these  curves  were  computed  will  now  be  described.* 


*Some  of  this  material  was  first  presented  by  the  author  as  a  paper  at  the  URSI-IRE  meeting  of 
Commission  2  at  Washington,  D.C.,  on  May  6,  1959.  This  work  was  later  described  in  an  NRL 
Report  (13)  and  in  correspondence  in  IRE  Transactions  (14).  The  results  there  described  were 
manually  computed.  Later  they  were  programmed  for  computer  solution  and  extended  to  higher 
frequencies;  this  work  was  described  at  the  1963  International  Symposium  of  the  IEEE  Group  on 
Antennas  and  Propagation  at  Boulder,  Colorado.  The  description  of  these  computations  is  presented 
here  for  the  first  time  in  report  form. 
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Van  Vleck’s  equations  (15-17)  for  absorption  per  nautical  mile*  of  oxygen,  7lt  and 
water  vapor,  y  ,  were  written  in  the  following  form: 


53,5  p 

y  -  - 

1  T2  X2 


h  a. 


h  a. 


1A2  +  h2  a*  (2-  1A)2  +  h2  a2 


+ 


h  a. 


(2  +  1A)2  +  h2  a  2 


(108) 


and 


7  = 
2 


8 . 5  x  104  p 
T5/2  IQ  278/TX2 


h  a„ 


(0.7418-  1A)2  +  h2  a  2 


h  a„ 


(0.7418+  1A)2  +  h2  a  2 


6.  3  ph  a3 
k2  T 


(109) 


The  symbols  are  defined  as  follows: 

P  -  total  atmospheric  pressure,  millibars 
T  -  absolute  atmospheric  temperature,  degrees  Kelvin 
\  -  radio  wavelength,  centimeters 

<*!  -  line-breadth  constant!  for  oxygen  at  sea  level,  taken  to  be  0.02  cm-1 

a2  -  line-breadth  constant  for  water-vapor  at  sea  level,  applicable  to  the  k  =  1.348 
cm  line  (f  =  22.2  GHz),  taken  to  be  0.1  cm"1 

a3  -  an  effective  line-breadth  constant  for  the  residual  absorption  due  to  water- 
vapor  lines  in  the  infrared  region,  taken  to  be  0.27  cm"1 

h  -  a  function  expressing  the  dependence  of  the  line-breadths  on  altitude  (pressure 
and  temperature) 

P  -  water-vapor  density,  grams  per  cubic  meter 

The  function  h  has  been  assumed,  following  Van  Vleck,  to  be  proportional  to  pres¬ 
sure  and  inversely  proportional  to  the  square  root  of  the  temperature.  Therefore,  with 
inclusion  of  a  constant  factor  to  give  h  =  latT  =  293°K  and  p  =  1013  mb, 

h  =  1. 7  X  10"2  pT"1/2  .  (110) 


Experiments  of  Hill  and  Gordy  (18)  indicate  that  the  dependence  on  temperature  may 
actually  be  about  t"°-8S  instead  of  T"0-5.  However,  this  would  result  in  a  maximum  in¬ 
crease  of  y  of  about  10  percent  at  the  altitude  of  lowest  atmospheric  temperature,  in  the 


*One  international  nautical  mile  equals  1.852  kilometers  exactly. 

tFollowing  Van  Vleck’s  practice  the  line-breadth  constant  is  here  expressed  in  terms  of  the  wave- 
number  breadth,  which  he  denoted  (tv/c).  The  value  0.02  is  equivalent  to  a  frequency  width  of  600 
MHz.  This  is  the  half-width  of  the  absorption  line  at  the  half-maximum  point. 
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30,000  to  80,000  foot  region.  The  increase  at  10,000  feet  is  about  3  percent,  and  of 
course  there  is  no  change  at  zero  altitude  where  the  absorption  is  greatest.  Since  the 
contribution  to  integrated  path  loss  above  30,000  feet  is  a  relatively  small  percentage  of 
the  total,  the  net  effect  of  a  change  to  the  Hill-Gordy  line-breadth  temperature  depend¬ 
ence  would  be  small. 

The  values  of  y2  and  y2 ,  and  the  total  y  -  y1  +  y2 ,  were  computed  at  altitude  inter¬ 
vals  of  1000  feet  up  to  a  100,000-foot  altitude,  assuming  the  pressure  and  temperature 
values  of  the  ICAO  Standard  Atmosphere  (19)  (essentially  the  same  as  the  ARDC  model 
atmosphere).  A  water- vapor  density  of  7.5  grams  per  cubic  meter  at  zero  altitude  was 
assumed,  and  the  height  dependence  was  patterned  after  the  average  statistics  for  Wash¬ 
ington,  D.C.,  over  a  10-year  period  (for  which  the  zero-altitude  average  density  was  6.18 
grams  per  cubic  meter  (13)).  These  atmospheric  models  were  entered  into  the  computer 
as  a  table  of  values  at  1000-foot  altitude  intervals. 

The  total  path  loss  for  two-way  propagation  (as  for  monostatic  radar)  was  computed, 
for  each  of  the  specified  elevation  angles,  by  integrating  27  as  a  function  of  distance  from 
the  antenna  along  the  actual  ray  paths.  The  one-way  loss  (as  for  radio  communication, 
telemetering,  etc.)  would  be  half  the  two-way  decibel  values.  The  range-height-angle 
relationship  was  computed  on  the  basis  of  the  CRPL  Exponential  Reference  Atmosphere 
(20)  (a  model  of  the  atmospheric  refractivity  profile),  for  a  surface  refractivity  Ns  =  313. 

However,  the  values  in  the  range  45  to  75  GHz  are  shown  dotted  in  Fig.  21,  Part  1, 
because  the  formulas  used  for  calculating  the  absorption  are  not  valid  in  the  region  of  the 
oxygen  resonance  at  60  GHz.  The  "resonance"  is  actually  a  group  of  several  closely- 
spaced  lines,  and  the  formula  employed  is  the  "centroid  approximation"  as  given  by 
Van  Vleck  (15),  which  is  valid  only  to  within  about  15  GHz  of  the  60-GHz  region.  The 
approximate  behavior  in  the  60-GHz  region  is  shown  dotted.  The  exact  behavior  has 
been  calculated  by  Meeks  and  Lilley  (21)  (for  oxygen  in  the  absence  of  any  water- vapor 
absorption)  by  summing  the  absorption  due  to  each  of  the  individual  resonance  lines. 

As  discussed  in  some  detail  by  Straiton  and  Tolbert  (22),  questions  exist  concerning 
the  exact  validity  of  the  Van- Vleck- Weisskopf  (VVW)  pressure-broadening  formula  (17) 
in  the  far  wings  (far  from  the  resonances).  The  use  of  the  value  as  =  0.27,  not  justifiable 
theoretically,  is  an  attempt  to  account  empirically  for  observed  deviations  from  theory 
of  the  residual  contribution  of  the  infrared  lines  to  the  water- vapor  absorption,  but  the 
correction  thus  obtained  is  not  fully  satisfactory.  Similarly,  oxygen  absorption  "in  the 
far  wings"  may  be  incorrect,  as  computed  from  the  VVW  formula.  However,  this  is  the 
best  theory  available  at  the  present  time,  and  there  is  considerable  reason  to  believe 
that  the  results  are  not  grossly  incorrect. 

The  results  are  in  good  agreement,  for  the  frequencies  at  which  they  made  computa¬ 
tions,  with  the  curves  published  by  Bean  and  Abbott  (23)  and  by  Hogg  (24)  for  total  transit 
of  the  troposphere.  (Hogg's  work  was  incidental  to  computations  of  tropospheric  noise 
temperature,  as  discussed  in  Sec.  6  of  this  report.)  To  the  author's  knowledge,  curves 
comparable  to  those  of  Fig.  22,  Part  1,  for  the  attenuation  between  points  within  the  tro¬ 
posphere,  have  not  been  calculated  elsewhere. 

The  curves  thus  computed  for  a  standard  atmosphere  are  of  course  subject  to  some 
variation  in  the  actual  atmosphere  when  the  pressure,  temperature,  and  water-vapor  con¬ 
tent  differ  from  the  standard  values  assumed.  However,  the  ordinary  variation  of  pres¬ 
sure  and  temperature  will  produce  negligible  variation  of  the  oxygen  absorption.  Since 
this  is  the  predominant  absorption  at  the  below  about  1000  MHz,  the  results  will  not  vary 
much  below  this  frequency.  From  about  3000  MHz  up  to  about  25  GHz,  however,  the 
water-vapor  absorption  is  a  significant  part  of  the  total  absorption,  and  consequently 
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variations  of  the  atmospheric  water-vapor  content  will  cause  significant  variations  of 
the  absorption.  If  it  is  desired  to  correct  for  this  variation,  it  is  necessary  to  plot  sep¬ 
arately  the  oxygen  and  water-vapor  absorptions,  instead  of  combining  them  as  was  done 
in  the  curves  of  Part  1.  The  two  absorption  components  are  directly  additive.  If  the 
water-vapor  component  is  plotted  for  a  standard  water-vapor  density  (e.g.,  7.5  g/m3), 
the  decibel  absorption  for  other  values  of  vapor  density  is  simply  proportional  to  the 
density.  The  extreme  range  of  water-vapor  density  is  from  about  2  g/m3  (in  the  arctic) 
to  about  20  g/m3  (in  the  tropics). 


Ionospheric  Absorption 

In  Part  1  the  absorption  that  occurs  in  the  ionosphere  was  discussed  briefly,  with 
reference  to  a  paper  by  Millman  (25).  Some  more  recent  work  by  Millman  (26)  provides 
additional  information.  In  Part  1  the  equation 

Li(dB)  =  Vf2  (1U) 

was  given  for  the  ionospheric  loss  in  decibels,  and  values  of  A  were  given  for  low-angle 
rays,  high-angle  rays,  daytime,  and  nighttime,  for  frequency  f  in  megahertz.  Although 
not  there  specified,  the  low-angle  figures  were  for  rays  of  10-degree  initial  elevation 
angle  and  the  high-angle  figures  were  for  90-degree  initial  angle.  Eq.  (Ill)  is  applicable 
to  radar  frequencies  at  and  above  100  MHz,  which  includes  the  great  majority  of  radars. 
The  attenuation  at  100  MHz  is  not  great,  and  because  of  the  inverse- square  frequency 
relationship,  ionospheric  absorption  is  unimportant  at  frequencies  higher  than  a  few 
hundred  megahertz. 

Millman's  more  recent  report  contains  curves  (Fig.  35)  from  which  values  of  A, 

Eq.  (Ill),  can  be  deduced  as  given  in  Table  9.  These  values  (for  10-  and  90-degree 
angles)  are  slightly  different  from  those  of  Millman's  earlier  paper.  It  must  be  realized, 
however,  that  the  values  depend  heavily  on  what  characteristics  of  the  ionosphere  are  as¬ 
sumed  to  be  "standard."  In  addition  to  the  pronounced  diurnal  variation  indicated  in 
Table  9,  there  are  also  lesser  seasonal  and  solar- sunspot- cycle  variations. 


Table  9 

Values  of  A  for  Eq.  (Ill) 


Elevation  Angle 
(degrees) 

A  (dB) 

Daytime 

Nighttime 

0 

2.1  x  104 

6.4  x  102 

5 

2.0  x  104 

6.0  x  102 

10 

1.8  x  104 

5.4  x  lo2 

30 

8.4  x  103 

3.0  x  io2 

90 

4.4  x  103 

1.6  x  102 

Absorption  by  Rain 

For  basic  radar  range  calculation,  no  rainfall  is  assumed,  since  in  most  areas  of 
the  world  rainfall  is  the  exception  rather  than  the  rule.  Nevertheless,  rain  does  occur 
often  enough  to  be  a  factor  in  assessing  overall  radar  detection- range  capability  at 
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microwave  frequencies.  Absorption  by  rain  is  seldom  important  below  about  10  GHz, 
however. 

The  computation  of  radar  maximum  range  for  a  specified  rainfall  attenuation  per 
unit  distance  and  a  specified  range  extent  of  the  rainfall  area  has  been  discussed  in  Part  1 
(Sec.  7).  The  problem  faced  by  a  radar  designer,  or  specification  writer,  is  to  specify 
the  characteristics  of  a  standard  rainfall  in  which  the  radar  must  meet  a  specified 
detection- range  capability.  This  is  a  problem  that  would  best  be  handled  by  a  meteorolo¬ 
gist;  but  to  the  author's  knowledge,  no  standards  of  this  type  have  been  proposed.  There¬ 
fore,  the  following  proposals  are  made  for  use  by  radar  designers  and  specification 
writers  until  more  authoritative  standards  become  available. 

A  light  rain  (27)  is  considered  to  be  one  of  about  1  mm/hr  (0.04  in/hr).  Such  rains 
can  extend  over  wide  geographical  areas,  and  to  altitudes  of  about  10,000  ft  (23).  Although 
the  rainfall  rate  may  be  anything  but  uniform  within  the  geographical  and  altitude  extent 
of  the  rain,  for  the  purpose  of  defining  a  convention  (see  Part  1,  Sec.  1)  for  absorption  by 
rain  it  is  not  unreasonable  to  assume  a  uniform  rate.  Therefore  it  is  proposed  that  a 
"standard  light  rain"  (for  computing  radar  performance)  be  defined  as  one  with  a  rate  of 
1  mm/hr  (0.04  in/hr)  extending  over  an  indefinitely  large  geographical  area  and  to 
10 ,000- ft  altitude. 

A  heavy  rain  (27)  is  one  of  about  16  mm/hr  (0.64  in/hr).  However,  such  rains  are 
usually  limited  in  their  geographical  extent,  seldom  exceeding  a  diameter  of  about  20 
miles.  Their  altitude  extent,  however,  may  be  20,000  ft  or  greater.  A  proposed  standard 
heavy  rain  is  therefore  one  of  16  mm/hr  covering  a  20- nautical- mile  diameter  geograph¬ 
ically,  and  extending  to  20 ,000- ft  altitude. 

An  intense  rain,  of  the  type  observed  in  severe  thunderstorms,  is  one  of  about  50 
mm/hr  (2  in/hr)  or  greater.  Rates  of  up  to  100  mm/hr  or  more  are  not  uncommon,  but 
these  are  usually  very  local  and  of  short  duration  in  nontropical  regions.  The  altitude  of 
a  thunderstorm  may  be  40,000  ft  or  greater.  A  standard  intense  rain  will  therefore  be 
arbitrarily  defined  as  one  of  50-mm/hr  rate,  2-nautical-mile  diameter,  and  40,000-ft- 
altitude  extent. 

It  is  to  be  emphasized  that  these  proposed  conventions  by  no  means  represent  or 
even  bracket  all  possible  rainfall  conditions.  Possibly  the  most  characteristic  feature 
of  rain  is  its  great  variability,  both  in  space  and  time.  The  assumption  of  uniform  rain 
rate  within  an  appreciable  volume  of  space  is  physically  unrealistic;  but,  for  computing 
absorption  such  an  assumption  is  justified  since  the  total  path  attenuation  is  proportional 
to  the  average  absorption  over  the  path. 

The  attenuation  produced  by  rain  of  specified  intensity  is  given  theoretically  by  J.W. 
Ryde  (28,29).  A  comparison  of  Ryde's  theory  with  numerious  experimental  results  has 
been  made  by  Medhurst  (30).  (Medhurst  also  applied  machine  computation  to  Ryde's 
equations,  thus  refining  the  basis  for  comparison  of  theory  and  experiment.)  He  found 
that  the  agreement  is  not  quantitatively  good,  but  there  is  qualitative  agreement.  The  ex¬ 
perimental  attenuations  are  often  greater  than  those  predicted  theoretically.  This  may  be 
partly  explained  by  inadequacies  of  the  experimental  techniques.  Medhurst  suggests  a 
number  of  possible  theoretical  factors  not  considered  by  Ryde's  theory  which  may  also 
partially  explain  the  discrepancies.  Semplak  and  Turrin  have  made  an  additional  sug¬ 
gestion  (31). 

Despite  this  disagreement,  the  Ryde  theory  is  probably  adequate  for  computing  ex¬ 
pected  attenuation  for  the  rainfall  conventions  proposed  here,  since  they  are  only  very 
approximately  related  to  reality  in  any  specific  case.  Medhurst  has  computed  values  of 
the  theoretical  attenuation  for  a  number  of  frequencies  as  a  function  of  the  precipitation 
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rate.  His  results  are  plotted  in  Fig.  4  for  the  rainfall  rates  of  the  proposed  conventions 
and  a  temperature  of  20 °C.  Temperature  affects  the  attenuation  in  a  very  complicated 
way.  Medhurst  (30)  has  published  a  table  of  temperature  correction  factors  based  on  the 
theory  of  Ryde. 


2000  10000  100000 

FREQUENCY*  MHZ 

Fig.  4  -  Attenuation  by  rain  of  specified  rates, 
after  Medhurst  (30) 


An  excellent  discussion  of  rain  attenuation  (and  also  of  radar  echoes  from  rain)  is 
contained  in  Ref.  1,  Chapters  7  and  8.  Although  it  was  written  three  decades  ago,  this 
book  is  still  a  most  valuable  source  of  information  on  these  subjects.  Absorption  and 
radar  reflection  by  hail,  snow,  and  fog  or  clouds  are  also  discussed.  Much  useful  mete¬ 
orological  information  and  radar  theory  are  also  given  by  Battan  (32). 


6.  NOISE  TEMPERATURE  CALCULATION 
System  Noise  Temperature 

The  objective  of  noise  temperature  calculation  is  to  obtain  a  measure  of  the  noise 
power  level  at  the  output  of  the  system.  This  in  turn  permits  computation  of  the  output 
signal-to-noise  ratio,  which  is  the  principal  factor  that  determines  the  system  sensitiv¬ 
ity.  It  is  not,  however,  the  sole  factor  if  the  term  signal-to-noise -ratio  is  interpreted  to 
mean  the  "instantaneous"  (or  single-pulse)  signal-to-noise  ratio.  The  ultimate  sensitivity 
is  then  also  affected  by  signal-processing  procedures  such  as  integration  or  averaging. 
(These  two  terms  are  virtually  synonymous.)  Here  the  term  signal-to-noise  ratio  will  be 
meant  in  the  instantaneous  or  preprocessing  sense. 

The  concept  of  system  noise  temperature  and  its  relationship  to  the  noise  tempera¬ 
tures  of  the  system  components  have  been  treated  in  some  detail  in  Part  1,  for  the  case 
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of  a  "single  response"  system;  multiple- response  systems  present  special  problems 
which  are  discussed  in  Refs.  33  and  34.* 

Most  radar  receiving  systems  are  of  the  single- response  type,  although  occasionally, 
at  the  very  highest  microwave  frequencies,  a  superheterodyne  receiver  without  preselec¬ 
tion  may  be  used.  Such  a  receiver  has  two  input  response  frequencies  —  the  desired  signal 
frequency,  and  the  image  frequency  (the  difference  and  the  sum  of  the  local  oscillator  and 
intermediate  frequencies).  For  the  appropriate  method  of  calculating  the  system  noise 
temperature  in  this  case,  Ref.  34  should  be  consulted;  additional  details  will  be  found  in 
Ref.  33. 


The  definitions  of  Ref.  34  do  not  contain  the  term  "system  noise  temperature."  It  is 
used  here  because  it  has  become  established  by  wide  usage  among  radio  and  radar  engi¬ 
neers;  it  is  a  de  facto  standard  term.  The  equivalent  term  used  in  Ref.  34  is  "operating 
noise  temperature  (of  a  system),"  designated  Top. 


If  a  receiving  system  consists  of  a  cascade  of  twoport  transducers,  it  is  in  principle 
possible  to  dissect  it  at  any  junction  of  two  of  the  system  elements,  and  then  to  regard 
the  portion  at  the  input  end  as  a  "source,"  and  the  remaining  portion  as  "the  receiver." 
The  "system"  then  has  two  elements,  and  the  following  equation  of  Ref.  34  applies: 


T  =  T;  +  T 

op  1  e 


(112) 


where  Tj  is  the  noise  temperature  of  the  source  (called  the  "input  termination"  in  Ref. 
34),  and  Te  is  the  effective  input  noise  temperature  of  the  receiver.  The  latter  quantity 
is  related  to  the  receiver  noise  factor  Fn  by  the  equation 


Te  =  (Fn  -  1)  T0  • 


(113) 


where  by  IEEE  standards  T0  =  290  °K. 


In  this  type  of  procedure,  the  dissection  is  essentially  arbitrary.  For  example,  if 
the  actual  system  consists  of  three  cascaded  components,  as  shown  in  Fig.  5,  it  could  be 
dissected  either  at  the  receiver  terminals  J2  (Case  1)  or  at  the  antenna  terminals  jx 
(Case  2).  In  the  former  case,  the  transmission  line  (which  generates  noise  because  of 
its  loss  factor  Lr  and  nonzero  thermal  temperature  (Tt)  is  regarded  as  a  part  of  the  an¬ 
tenna,  and  the  intrinsic  antenna  noise  temperature  Ta  is  modified  by  the  effect  of  the  line 
loss.  In  the  latter  case,  the  receiver  noise  temperature  is  modified  by  the  line  loss.  In 
either  case,  the  output  signal-to-noise  ratio  can  be  correctly  computed,  but  the  value  of 
Top  obtained  in  the  two  cases  will  be  different.  Since  this  point  sometimes  causes  con¬ 
siderable  confusion,  it  will  be  discussed  here  in  some  detail.  The  notation  Top  (operat¬ 
ing  noise  temperature)  will,  however,  now  be  replaced  by  the  more  commonly  used  Ts 
(system  noise  temperature). 


The  noise  temperature  Tj  in  Eq.  (112)  multiplied  by  Boltzmann's  constant  k  =  1.38  x 
10'23  (watt- sec/ °K)  and  by  the  noise  bandwidth  of  the  system  Bn  gives  the  available  out¬ 
put  noise  power  of  the  input  termination  —  i.e.,  the  noise  power  that  would  be  delivered 
to  a  matched  load  of  bandwidth  Bn  connected  at  the  output  terminals  of  the  source/  When 


*The  problem  of  multiple-response  systems  is  considered  in  appreciably  greater  detail  in  Ref.  33 
than  in  the  IEEE  Standard  and  tutorial  paper,  Ref.  34.  However,  the  definitions  used  in  Ref.  33, 
which  was  written  before  the  IEEE  Standard  was  issued,  do  not  conform  exactly  to  the  IEEE  Stand¬ 
ard  definitions.  Since  definitions  are  basically  arbitrary,  the  treatment  in  Ref.  33  is  valid  and  use¬ 
ful  in  terms  of  its  own  definitions,  but  it  differs  in  some  details  from  the  treatment  of  Ref.  34. 
f This  formula  is  actually  an  approximation,  analogous  to  the  Rayleigh-Jeans  approximation  of 
Planck's  radiation  law.  The  more  exact  expression  is  hfBn/(  eh f  /k/  -  1),  where  h  is  Planck's 
constant  and  f  is  frequency,  Hz.  The  approximation  is  invalid  if  f/T  >  4  x  10 8. 
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Fig.  5  -  A  three-component  cascade  receiving  system 


the  system  is  dissected  at  the  antenna  terminals  (Case  2)  it  is  evident  that  the  correct 
value  for  Tj  is  thus  simply  Ta,  the  antenna  noise  temperature. 

If  the  dissection  is  at  the  receiver  terminals  (Case  1),  however,  the  value  of  Tj  re¬ 
flects  the  resultant  available  noise  power  (at  the  receiver  terminals)  of  both  the  antenna 
and  the  transmission  line.  As  has  been  discussed  in  Part  1,  the  available  output  noise 
power  in  bandwidth  Bn  of  a  passive  twoport  transducer  (e.g.,  transmission  line)  of  avail¬ 
able  loss  factor  Lr  and  thermal  temperature  Tt  is 

P„  =  kTt(l-  1/Lr)  Bn  .  (114) 

This  figure  assumes  that  an  input  termination  (e.g.,  antenna)  is  connected  to  the 
transmission- line  input  terminals,  but  it  does  not  include  the  noise  power  at  the  receiver 
terminals  due  to  the  antenna.  This  contribution  is  directly  additive  to  the  transmission¬ 
line  noise  power,  as  mentioned  in  Part  1,  because  the  two  noise  sources  are  not  coherent 
(phase  correlated).  The  antenna  noise  power  at  this  point  is 

Po(ant)  =  kT8Bn/Lr  .  (115) 

The  divisor  Lr  expresses  the  fact  that  the  antenna  noise  is  attenuated  by  the  transmission 
line.  The  total  noise  is  therefore 


kT!Bn  =  kBn[Ta/Lr  +  Tt(l-  1/Lr>] 


(116) 


The  term  on  the  right  in  brackets  is  evidently  equal  to  T(.  (The  "prime"  on  Tj  is  here 
used  to  express  the  fact  that  it  is  not  the  same  as  Tj  when  the  system  is  dissected  at  the 
antenna  terminals.)  The  total  system  noise  temperature  for  Case  1  then  becomes 

Ts(l)  "  V^r  +  Tt(l-1/Lr)  +  Te  ■  (117) 

(by  substituting  the  solution  for  Tf  in  Eq.  (116)  for  Tj  in  Eq.  (112).) 

For  Case  2,  Tj  is  simply  Ta,  but  Te  (which  will  now  be  denoted  T')  must  express 
the  combined  effect  of  the  transmission  line  and  the  receiver.  The  principle  invoked  to 
solve  this  problem  is  that  the  product  kTaBnG  must  represent  the  contribution  of  the 
transmission-line-and-receiver  combination  to  the  total  system  output  noise  power, 
where  G  is  the  overall  power  gain  of  the  combination.  It  is  evident  that  this  overall  gain 
is  Gr/Lr  where  Gr  is  the  receiver  gain  and  Lr  is  the  transmission- line  loss  factor. 
Therefore  the  equation  that  must  hold  is 


kTeBnGr/Lr  -  PD(line)  +  Po  (  receiver) 


(118) 


where  PD  denotes  noise  power  at  the  receiver  output  terminals.  From  Eq.  (114)  it  is 
evident  that 
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Po(Hne)  =  kTt(l“  1/Lr)  BnGr 

and  it  is  obvious  that 

Po( receiver)  “  k  Te  Bn  Gr  ' 

By  substituting  Eqs.  (119)  and  (120)  into  Eq.  (118)  it  is  found  that 

t;  =  Tt(Lr-  1)  +  LrTe  . 


(119) 


(120) 


(121) 


Hence  the  total  system  noise  temperature  for  Case  2  is 


Ts(2)  =  Ta  +  Tt(Lr-l)  +  LrTe  • 
Comparing  this  result  with  Eq.  (117),  it  is  readily  seen  that 


(122) 


Ts(2)  -  Lr  Ts(l)  ' 


(123) 


It  will  now  be  shown  that  either  Eq.  (117)  or  (123)  can  be  used  to  calculate  the  system 
signal- to- noise  ratio  correctly.  Suppose  that  a  signal  power  Ps(ant)  is  available  at  the 
output  terminals  of  the  antenna.  The  signal  power  at  the  system  output  will  be 


Ps(out)  Ps(ant)  ' 

(124) 

In  Case  1  the  noise  power  output  is  given  by 

Pn  (out)  ”  k  Ts(l)  BnGr 

(125) 

and  in  Case  2  it  is  given  by 

Pn(out)  =  k  Ts  (  2  )  Bn  Gr'/Lr  ’ 

(126) 

But  since  Ts(2)  =  LrTs(1)  by  Eq.  (123),  these  two  values  of  Pn,out)  are  identical,  and  the 
signal- to- noise  ratio  PS(out)/pn(out)  is  therefore  the  same  in  both  cases,  as  was  to  be 
proved. 


Thus  it  makes  no  difference  which  dissection  point  is  chosen,  except  for  the  following 
consideration.  In  calculating  radar  maximum  range,  the  system  loss  factor  L  must  con¬ 
tain  as  one  of  its  factors  the  receiving  line  loss  factor  Lr  if  the  Case  1  representation  of 
Ts  is  used.  If  the  Case  2  representation  is  used,  however,  in  effect  there  is  no  receiving 
line  loss  (i.e.,  Lr  =  1),  because  the  receiving  line  is  considered  to  be  a  part  of  the  re¬ 
ceiver;  its  noise  is  incorporated  into  the  computation  of  the  effective  noise  temperature 
of  the  receiver,  and  its  loss  factor  Lr  reduces  the  effective  receiver  gain. 

The  Case  1  representation  was  employed  in  the  first  edition  (1962)  of  Part  1  of  this 
report,  and  Lr  was  included  in  the  computation  of  L.  In  the  second  edition  (1969),  the 
Case  2  representation  is  used,  and  the  factor  Lr  is  not  included  in  computation  of  L.  As 
the  foregoing  discussion  implies,  either  method  results  in  the  same  calculated  range. 

The  fact  is,  of  course,  that  in  a  system  consisting  of  a  multiplicity  of  cascaded  units, 
the  dissection  of  the  system  can  be  made  at  any  junction  of  two  system  components,  and 
in  this  way  many  possible  values  of  the  system  noise  temperature  can  be  obtained,  all 
valid,  and  all  applicable  to  one  and  the  same  system.  In  each  such  case,  however,  the 
receiving  line  loss  factor  Lr  must  be  evaluated  differently;  it  must  represent  only  those 
losses  that  occur  prior  to  the  "dissection"  point. 
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This  matter  is  discussed  in  some  detail  in  Ref.  33,  in  terms  of  what  is  there  called 
the  referral  principle.  The  system  noise  temperature  is  said  to  be  referred  to  the  point 
in  the  system  which,  in  the  foregoing  discussion,  is  called  the  dissection  point.  Very 
general  equations  are  there  presented  for  the  system  noise  temperature  —  i.e.,  for  the 
temperature  of  a  system  of  n  cascaded  components  with  the  system  noise  temperature 
referred  to  the  junction  of  the  mth  and  (m  +  l)th  components  of  the  system  (m  <  n). 

In  1962  there  seemed  to  be  historical  precedents  for  taking  the  receiver  input  ter¬ 
minals  as  the  system-noise-temperature  reference  point.  However,  in  the  1962  edition 
of  Part  1  of  this  report,  the  concept  of  the  system  input  noise  temperature  was  also  dis¬ 
cussed.  At  that  time,  the  author  also  chose  to  define  antenna  noise  temperature  in  terms 
of  the  noise  power  captured  from  space  by  the  antenna  aperture  prior  to  its  reduction  by 
any  ohmic  losses  in  the  antenna;  this  decision  was  made  in  the  absence  of  any  IEEE 
standard  definition  of  antenna  noise  temperature.  Subsequently  (35)  an  IEEE  standard 
definition  has  been  adopted;  it  defines  the  antenna  noise  temperature  in  terms  of  the 
noise  power  available  at  the  actual  terminals  of  the  antenna.  As  discussed  in  Part  1, 
this  noise  is  a  composite  of  noise  received  from  external  radiating  sources  and  noise 
generated  in  dissipative  portions  of  the  antenna. 

Since  1962  there  has  been  a  trend  toward  defining  the  system  noise  temperature  with 
the  antenna  terminals  as  the  reference  point.  Therefore  this  practice  has  been  followed 
in  the  second  edition  of  Part  1,  in  terms  of  the  IEEE  definition  of  the  antenna  noise  tem¬ 
perature. 

For  the  special  case  of  this  reference  point,  the  general  equation  for  the  system 
noise  temperature  of  a  single- response  cascade  system  of  N  twoport  transducers  (not 
including  the  antenna,  which  is  regarded  as  the  input  termination  of  the  first  transducer) 
becomes 


Ts  -  Ta  +  £  Te(i)/Gi  .  (127) 

i=  1 

which  was  given  in  Part  1  as  Eq.  (35).  In  this  equation  Te(i)  is  the  effective  input  noise 
temperature  of  the  ith  transducer  and  cq  is  the  available  gain  of  the  system  between  its 
input  terminals  (antenna  terminals)  and  the  input  terminals  of  the  ith  transducer.  As 
noted  in  Part  1,  by  this  definition  Gl  =  l  always.  For  the  three- component  (two- 
transducer)  system  that  has  been  analyzed  here,  G2  =  l/Lr,  Te(1)  =  Tt(Lr-i),  and 
Te(  2>  =  Te ,  where  Te  without  additional  subscripts  denotes  the  receiver  noise  tempera¬ 
ture  as  defined  by  Eq.  (113)).  Thus  Eq.  (127)  becomes  Eq.  (122).  For  further  details 
and  explanations,  see  Ref.  33. 


Antenna  Noise  Temperature 

As  discussed  under  the  preceding  heading,  the  IEEE  standard  definition  of  antenna 
noise  temperature  defines  it  to  be  the  composite  result  of  noise  from  external  radiating 
sources  and  noise  from  the  internal  dissipative  parts  of  the  antenna.  In  Part  1,  curves 
are  presented  (Fig.  11)  which  give  the  noise  temperature  of  a  typical  radar  antenna  due 
to  external  radiating  sources,  for  conventional  assumptions  concerning  the  external  en¬ 
vironment,  and  omitting  the  contribution  of  radiation  from  the  ground  (earth).  The  text 
of  Part  1  explains  how  the  value  of  Ta  obtained  from  these  curves  is  modified  to  take 
into  account  the  ground  radiation  and  the  effect  of  antenna  dissipative  losses.  The  pur¬ 
pose  of  this  discussion  will  be  to  explain  in  some  detail  how  the  curves  of  Fig.  11,  Part  1, 
were  calculated. 
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As  there  explained,  these  curves  represent  the  composite  effect  of  (i)  galactic  noise, 
(ii)  sun  noise,  (iii)  tropospheric  noise,  and  (iv)  cosmic  blackbody  radiation.  These  con¬ 
tributing  sources  are  in  general  functions  of  the  radar  frequency  and  the  antenna  beam 
pointing  direction.  To  plot  curves  that  can  be  used  as  conventions  for  radar  range  cal¬ 
culation,  the  variation  of  galactic  noise  with  galactic  coordinates  is  ignored,  and  a  mean 
value  is  assumed  (actually  the  geometric  mean  of  maximum  and  minimum  values).  The 
variation  of  tropospheric  noise  with  antenna  beam  elevation  angle  is  taken  into  account, 
so  that  a  family  of  curves  with  beam  elevation  angle  as  the  family  parameter  results. 

The  sun  is  assumed  to  be  in  the  visible  sky,  but  not  in  the  main  beam  of  the  antenna  pat¬ 
tern;  a  conventional  value  of  sidelobe  level  in  the  sun  direction  is  assumed  (unity-gain  or 
0-dB  side  lobes).  The  cosmic  blackbody  radiation  is  assumed  to  be  independent  of  both 
frequency  and  direction  and  to  be  of  noise  temperature  2.7 °K  (36).  This  small  contribu¬ 
tion  to  the  total  antenna  temperature  is  of  little  importance  ordinarily,  but  it  may  be  of 
significance  for  antennas  of  the  ultralow- noise  variety  (37).  (The  value  3.5 °K  was  orig¬ 
inally  estimated  for  this  temperature,  but  the  value  most  recently  reported  is  2.7°K.) 

The  well-known  basic  equation  for  the  noise  temperature  of  an  antenna  is 

277  77 

Ta  =  —  f  j  T(0,0)  G(0,<£)  sini9d0d0  ,  (128) 

47/  J0 

in  which  T (6,4>)  is  the  noise  temperature  "seen"  by  the  antenna  in  the  direction 
and  G (9,<p)  is  the  antenna  gain  in  that  direction.  Since  T(0,<£)  is  seldom  known  in  suffi¬ 
cient  detail  to  permit  this  integration  to  be  carried  out  exactly,  it  is  usual  to  average 
over  finite  segments  of  solid  angle.  If  the  entire  477  steradians  is  thus  sub¬ 
divided  into  n  segments,  Eq.  (128)  becomes 


=  L 

i  -  1 


Aa(i) 


(129) 


where  Ta(i)  is  the  average  noise  temperature  within  the  ith  solid-angle  segment,  whose 
solid-angular  size  will  be  denoted  fi4.  A  general  formula  for  Ta(i^  is 


fii  G.  t. 


>(i) 


47 T 


-  a-  T-  , 


(130) 


where  Gi  and  T*  are  the  average  values  of  G(0,<£)  and  T(@,<£)  within  fii.  The  fraction 
oii  =  fii  Gi/4 77  can  be  interpreted,  in  terms  of  a  transmitting  antenna,  as  the  fraction  of 
the  total  power  radiated  into  the  solid  angle  fii . 

At  some  frequencies  the  noise  temperature  in  a  particular  direction  is  partially  due 
to  astronomical  sources  (the  galaxy  and  the  sun)  and  partly  due  to  radiation  by  the  lossy 
atmosphere.  In  such  cases  the  net  temperature  is  the  noise  temperature  of  the  distant 
source  reduced  by  the  loss  factor  of  the  atmosphere  and  augmented  by  the  atmospheric 
radiation.  For  example,  in  the  case  of  the  galactic  noise  the  formula  is 


T 


a(gal) 


gal 


(131) 


where  Ta(gal)  is  the  contribution  of  the  galactic  solid  angle  to  the  total  antenna  tempera¬ 
ture,  agal  is  defined  by  Eq.  (130),  Tgal  is  the  average  galactic  noise  temperature  before 
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attenuation  by  the  atmosphere,  and  Tatm  is  the  atmospheric  noise  temperature.  A  similar 
procedure  is  followed  for  the  sun  noise  contribution. 

As  explained  above,  the  ground  noise  contribution  was  assumed  to  be  zero  in  plotting 
Fig.  11  of  Part  1.  The  procedure  there  described  for  adding  a  ground  noise  contribution 
simply  corresponds  to  adding  a  term  aG  Tc  in  Eq.  (129),  where  aG  represents  the  fraction 
of  4?r  steradians  subtended  by  the  ground  as  viewed  by  the  antenna  multiplied  by  the  av¬ 
erage  gain  of  the  side  lobes  that  view  the  ground,  and  TG  is  the  noise  temperature  of  the 
ground. 

If  the  ground  is  completely  absorbing,  its  noise  temperature  is  equal  to  its  thermal 
temperature,  TtG  (Kelvin).  If  it  is  completely  reflecting,  its  noise  temperature  is  deter¬ 
mined  by  its  reflection  of  the  galactic  noise,  sun  noise,  etc.  In  general,  its  emissivity  is 
6=1-$,  where  £R  is  its  reflectivity,*  or  the  fraction  of  incident  power  that  is  reflected. 
In  these  terms, 


TG  =  STtG.  (132) 

The  galactic  noise  is  assumed  to  follow  the  frequency  law  f~ 2 • 5  (38,  39),  and  its 
noise  temperature  at  100  MHz  is  taken  to  be  3050  °K.  This  is  the  geometric  mean  of  the 
maximum  and  minimum  values  reported  by  Brown  and  Hazard  (40),  which  at  100  MHz  are 
18,650  °K  and  500 °K.  The  geometric  mean  was  taken  so  as  not  to  give  undue  weight  to 
galactic  "hot  spots”  of  small  angular  dimensions.  Also,  this  procedure  causes  the  mean 
to  fall  halfway  between  the  maximum  and  minimum  values  (dashed  lines)  on  the  logarithmic 
temperature  scale  of  Fig.  11,  Part  1.  The  formula  thus  obtained  for  the  galactic  noise 
contribution  to  the  antenna  temperature  is 

Ta(gal)  =  “gal  T100  ( 100/f MHz) 2 ' S/Latm  >  (133) 

where  T100  is  the  galactic  temperature  at  100  MHz,  fMHz  is  the  frequency  in  megahertz, 
ff gai  is  given  by  Eq.  (130),  and  Latm  is  the  atmospheric  loss  factor  (a  function  of  eleva¬ 
tion  angle) . 

The  value  assumed  for  ctgal  in  computing  the  curves  of  Fig.  11,  Part  1,  is  unity, 
because  these  curves  are  for  an  ideal  antenna  which  has  no  side  or  back  lobes  that  view 
the  ground.  Correction  factors  to  be  applied  for  actual  antennas  are  discussed  in  Part  1. 

A  typical  value  of  agal  for  a  real  antenna  is  about  0.85. 

The  noise  temperature  of  the  quiet  sun,  Tsun,  was  taken  from  a  curve  published  by 
Kuiper  (41),  after  Smerd  (42).  The  values  in  the  range  of  frequencies  plotted  in  Fig.  11, 
Part  1,  are  given  in  Fig.  6.  The  sun  noise  contribution  to  the  antenna  temperature  was 
then  calculated  from  Eq.  (130)  with  asun  (=  fisun  G/4tt)  equal  to  4.75  X  10"6,  with  G  =  1. 

The  sun’s  visible  (optical)  diameter  subtends  an  angle  of  1/2  degree  as_seen  from  the 
earth,  and  this  value  of  «sun  corresponds  to  this  diameter.  The  value  G  =  1  means  that 
the  sun  was  assumed  to  be  viewed  in  a  unity-gain  (0-dB)  side  lobe  of  the  antenna  pattern. 

It  was  also  assumed  that  the  sun  was  in  an  active  state,  so  that  its  noise  temperature  was 
ten  times  the  quiet  level.  Consequently,  the  formula  for  the  sun  noise  contribution  to 
total  antenna  temperature  is 


*It  is  to  be  noted  that  the  reflectivity  31  is  not  the  same  as  the  specular  reflection  coefficient  p  dis¬ 
cussed  in  Part  1,  Sec.  6.  Where  p  is  the  ratio  of  the  specularly  reflected  field  strength  (e.g.,  volts 
per  meter)  to  the  incident  field  strength,  31  is  the  ratio  of  the  total  power  (e.g.,  watts)  reflected  from 
a  unit  area  of  surface,  including  both  specular  and  diffuse  reflection,  to  the  power  incident  on  the 
area. 
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Fig.  6  -  Calculated  apparent  noise 
temperature  of  the  quiet  sun,  after 
Smerd  (42).  A  coronal  temperature 
of  10 6  degrees  and  a  chromospheric 
temperature  of  104  degrees  Kelvin 
were  assumed.  Dots  represent  ob¬ 
servational  data  taken  during  the 
period  1942-1950. 


FREQUENCY  (MHl) 


ra(.Un)  =  4.75  X  10-5  Tsun/Latn 


Although  the  sun's  optical  diameter  is  1/2  degree,  its  "noise  diameter"  is  consid¬ 
erably  greater  at  the  lower  end  of  the  radar  frequency  spectrum,  because  the  corona  is 
there  the  principal  source  of  radiation, 'while  at  the  upper  end  of  the  spectrum  the  radia¬ 
tion  comes  from  the  chromosphere.  This  increase  of  effective  diameter  at  the  lower 
frequencies  is  taken  into  account  by  the  method  of  calculating  Tsun  (41). 


The  atmospheric  loss  factors  are  obtainable  from  Fig.  21  of  Part  1.  Since  the 
values  shown  there  are  decibel  values  for  two-way  transit  of  the  troposphere,  the  rela¬ 
tion  between  the  atmospheric  loss  factor  Latm  in  Eqs.  (131),  (133),  and  (134)  and  the  dec¬ 
ibel  (dB)  values  of  Fig.  21,  Part  1,  is 


Latm  =  10dB/2°  • 


The  noise  temperature  of  the  atmosphere  Tatm  is  given  (43)  by  the  equation 


rm  -0.2303J  7(r)dr 

Tatm=  0.2303  J  ?Tt(atra)  e 


in  which  the  variable  of  integration  R  is  distance  along  the  ray  path  from  the  antenna, 
y  (R)  is  the  attenuation  in  decibels  per  unit  distance  along  the  path,  and  Tt(atm)  is  the 
thermal  temperature  of  the  atmosphere.  (The  numerical  factor  0.2303  occurs  in  this 
equation  because  7  is  defined  to  be  in  decibels  rather  than  nepers  per  unit  distance.) 

The  values  of  7  were  computed  as  described  in  Sec.  5,  and  the  Tt(atm)  values  used  were 
those  of  the  ICAO  Standard  Atmosphere  (19). 


These  results  were  applied  to  Eq.  (129)  to  compute  antenna  noise  temperature  as  a 
function  of  frequency  from  100  MHz  to  100  GHz,  for  seven  elevation  angles.  The  compu¬ 
tations  were  done  by  the  NRL  CDC-3800  computer  from  a  FORTRAN  program  (TAPLOT) 
written  by  the  author,  and  machine  plotted  by  the  NRL  Gerber  Model  875  Automatic  Draft¬ 
ing  Machine  from  punched  paper  tape  obtained  from  the  computer  results.  The  resulting 
plot  is  Fig.  11  of  Part  1.  In  the  left  part  of  the  figure  the  dashed  lines  represent  maxi¬ 
mum  and  minimum  values  of  the  galactic  and  solar  noise  temperatures,  and  in  the  right 
part  they  represent  maximum  (low- angle)  and  minimum  (high- angle)  tropospheric  noise. 
The  maximum  curve  assumes  a  sun  noise  temperature  100  times  the  quiet  level  of  Fig.  6, 
while  the  minimum  curve  assumes  zero  sun  noise  (nighttime  radar  operation).  The  slight 
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spreading  of  the  solid  curves  in  the  UHF  region  is  due  to  the  variation  of  absorption  of 
galactic  noise  by  atmospheric  loss  (see  Eq.  (131)),  causing  a  higher  antenna  temperature 
at  higher  elevation  angles.  That  is,  the  reduction  of  the  galactic  noise  at  the  low  angles 
caused  by  atmospheric  absorption  is  greater  than  the  attendant  tropospheric  noise  radia¬ 
tion,  in  this  frequency  region.  The  peculiar  bump  in  the  curves  at  about  500  MHz  is  due 
to  the  knee  of  the  sun  noise  curve  at  about  that  frequency  (Fig.  6). 


Transmission- Line  Noise  Temperature 
In  Part  1,  Sec.  5,  the  formula 


Pno  =  kTtBn(l-  1/L)  (137) 

was  given  for  the  noise  power  output  in  a  bandwidth  Bn  of  a  twoport  passive  transducer 
(e.g.,  transmission  line)  whose  thermal  temperature  is  Tt  and  whose  loss  factor  is  L. 
This  results  in  the  formula 


Tr=Tt(L-l)  (138) 

for  the  input  noise  temperature  of  the  transducer.  Eq.  (137)  was  derived  by  Dicke  (43) 
in  the  following  manner. 

The  derivation  will  be  made  in  terms  of  Fig.  5,  in  which  all  components  are  assumed 
to  be  linear  and  passive.  Initially  it  is  assumed  that  the  source  is  at  thermal  tempera¬ 
ture  Tt,  and  the  other  two  components  are  at  zero  temperature,  so  that  they  generate  no 
noise  power.  In  that  case,  the  available  noise  power  at  is  kTtBn  (in  accordance  with 
Nyquist's  Theorem),  and  that  at  J2  is  kTtBn/L,  where  L  is  the  loss  factor  (available  loss) 
of  the  twoport  transducer.  It  is  not  important  for  this  derivation  that  if  the  transducer 
were  at  zero  temperature  its  dissipative  loss  would  be  zero.  In  fact,  L  is  the  available 
loss,  which  is  the  product  of  the  mismatch  loss  at  Jj  and  the  dissipative  loss,  if  any,  in 
the  transducer.  But  the  argument  being  given  is  a  primarily  logical  rather  than  a  purely 
physical  one,  and  it  is  permissible  here  to  assume  that  the  loss  factor  is  independent  of 
the  temperature. 

If  now  the  temperature  of  the  twoport  transducer  is  raised  to  the  temperature  Tt, 
all  other  temperatures  remaining  unchanged,  the  available  power  at  J2  will  be  kTt  Bn, 
again  by  Nyquist's  Theorem.  Since  the  contribution  of  the  source  to  the  power  at  J2  must 
be  the  same  as  before,  the  contribution  of  the  twoport  transducer  must  be  the  difference 
between  the  available  power  at  J2  before  and  after  the  temperature  was  raised.  This 
difference  is 


Pno  =  kTtBn  -  kTtVL  =  kTtBn(l-  1/L) 


(139) 


as  was  to  be  proved. 

This  result  was  derived  by  assuming  that  the  source  and  the  transducer  are  brought 
to  the  same  thermal  temperature.  However,  because  the  system  is  linear  and  the  noise 
voltages  in  the  source  and  transducer  are  noncoherent,  it  is  evident  that  the  noise  con¬ 
tribution  of  the  transducer  to  the  total  noise  at  J2  is  not  dependent  on  the  temperature  of 
the  source.  Consequently  the  result  deduced  by  assuming  a  special  temperature  for  the 
source  is  actually  a  general  result. 

Further  details  of  this  matter  and  other  matters  discussed  in  this  section  will  be 
found  in  Ref.  33. 
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7.  DERIVATION  OF  FORMULA  FOR  RANGE  PREDICTION  ERROR 

In  Part  1,  Eq.  (92)  is  presented  without  proof.  It  expresses  the  range  prediction 
error,  or  standard  deviation,  as  a  function  of  the  standard  deviations  of  the  factors  in 
the  range  equation.  This  equation  was  derived  in  1961  by  John  Wood,  then  of  the  NRL 
Radar  Division  Mathematics  Staff,  in  response  to  a  request  by  the  author  to  S.  F.  George, 
Head  of  the  Mathematics  Staff.  Wood's  derivation  is  presented  here.  The  starting  point 
was  the  following  radar  range  equation,  which  was  given  in  Part  1: 


R  =  k 


pl/4Tl/4Gl/4Gl/4crl/4F 

f  1/2  T1/4v1/4t  1/4 

s  0 


(140) 


The  symbols  are  defined  as  follows: 

k  -  a  constant  depending  on  the  units  used 
Pt  -  transmitter  pulse  power 
r  -  pulse  length 
Gt  -  transmitting  antenna  gain 
Gr  -  receiving  antenna  gain 
a-  -  target  cross  section 
F  -  pattern- propagation  factor 
f  -  radar  frequency 
Ts  -  system  noise  temperature 

v0  -  visibility  factor  (minimum- detectable  signal- to- noise  power  ratio) 

L  -  system  loss  factor. 

The  problem  as  presented  was  to  find  the  standard  deviation  of  R,  given  the  individual 
standard  deviations  of  the  parameters  on  the  right-hand  side  of  Eq.  (139).  The  right-hand 
side  contains  the  product  (and  quotient)  of  terms.  Convenient  theorems  exist  concerning 
the  standard  deviation  of  the  sum  (and  difference)  of  terms.  The  suggestion  was  there¬ 
fore  made,  in  presenting  the  problem  to  Wood,  that  it  would  be  tractable  in  terms  of  a 
logarithmic  formulation.  This  approach  was  followed. 

The  following  theorems  are  employed.  Given  a  function  f  (xr  x2,  ....  xn),  it  is  known 

that 


(§  f)2  = 


(141) 


For  the  case  when 


f  (X1>  x2-  •  •  •  xn)  =  VlOl)  y2(x2>  •  •  ■  yn(xn) 


(142) 
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the  corresponding  logarithmic  relationship  is 


n 

1oe  f  =  2Z  log  yi(xi)  • 

i=  1 


Differentiating, 


3  log  f  _  _l  _3f  _  _l_  3yi 

3xj  f  y j  Bxj 


and  Eq.  (141)  therefore  gives 


<11) 1 

f2 


E 

i-  1 


<®xi)2 


(143) 


(144) 


(145) 


The  left-hand  side  of  this  equation  is  the  square  of  the  fractional  error  in  f.  Apply¬ 
ing  the  right-hand  side  of  Eq.  (144)  to  Eq.  (140)  gives 


I  11  =  _J_  I  p-3/4  -  _J_  . 
R  3Pt  pJ/4  4  *  4Pt 


(146) 


This  shows  that  an  expression  can  be  obtained  for  the  fractional  error  of  r  in  terms  of 
the  fractional  errors  of  the  other  parameters.  The  result  is 


(147) 


This  result  is  equivalent  to  Eq.  (92),  Part  1. 
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Appendix 


GRAPH  FOR  PREDICTING  THE  EFFECT  OF  A 
FINE -LOBE  INTERFERENCE  PATTERN 


The  graph  shown  below  is  a  reproduction  of  Fig.  1  of  Ref.  9,  which  is  discussed  in 
Sec.  4  of  this  report.  It  is  reproduced  here  because  Ref.  9  is  no  longer  readily  available. 
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R/R50  RATIO  OF  RANGE  TO  RANGE  OF  50%  BLIP-SCAN  RATIO  IN  FREE  SPACE 


Fig.  1  -  The  blip-scan  ratio,  averaged  over  one  lobe  of  the  interference  pattern,  as  a  function 
of  the  normalized  range.  Rso  ,  the  free-space  range  for  a  50-percent  blip-scan  ratio  on  a  given 
target,  is  equal  to  R0  v^r~ ,  where  a-s0  is  the  median  cross  section  of  the  target  and  rq  the 
free-space  range  on  a  one-square-meter  nonscintillating  target.  D  is  divergence  factor. 
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